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ABSTRACT 



Context. The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through 
the circumstellar disk onto the central star. 

Aims. We constrain the masses of the envelopes, disks, and central stars of a sample of low-mass protostars and compare the results 
to theoretical models for the evolution of young stellar objects through the early protostellar stages. 
Methods. A sample of 20 Class and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular 
resolution (typically 2") with the Submillimeter Array. Using detailed dust radiative transfer models of the interferometric data, as well 
as single-dish continuum observations, we have developed a framework for disentangling the continuum emission from the envelopes 
and disks, and from that estimated their masses. For the Class I sources in the sample HCO + 3-2 line emission has furthermore been 
observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, making it possible to determine the 
masses of the central stars. In the other sources the disks are masked by optically thick envelope and outflow emission. 
Results. Both Class and I protostars are surrounded by disks with typical masses of about 0.05 M , although significant scatter is 
seen in the derived disk masses for objects within both evolutionary stages. No evidence is found for a correlation between the disk 
mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from ~ 1 M Q in the 
Class stage to < 0.1 M e in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses 
in the Class stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars 
contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the 
embedded protostellar stages, with most of the material from the ambient envelope accreted onto the central star. Theoretical models 
. tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. 

Conclusions. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution 
(although these disks are not necessarily rotationally supported) and rapidly process material accreted from the larger scale envelope 
onto the central star. 

Key words, star: formation - circumstellar matter - planetary systems: protoplanetary disks - radiative transfer 



1 . Introduction (e. g.. lAndre et all 1 1 993l |2000|) . but few hard constraints on the 

actual evolution of the matter through the embedded protostel- 

The fundamental problem in studies of the evolution of low- ^ ^ Direct measurements of the steUar> diski and 

mass young stellar objects is how material is accreted from the envel masses duri these ivotal st can lace individ . 
larger scale envelope through the protoplanetary disk onto the ^ tostars m their evolutionary context by address- 
central star - and, in particular, what evolutionary time-scales . quest ions such as what fraction of the initial core mass has 
determine the formation and growth of the disk and star and been ontQ ±e Qemal what fraction ig bd carried 
the dissipation of envelopes. Our theoretical understanding of fe ^ action of ±e tostellar outflows and what frac . 
the evolution of YSOs has typically relied on the coupling be- don hag (sQ . far) ended ^ a difjk In Mf , we anaJ 
tween theoretical studies and empirical schemes based, for ex- si le _ dish and i nter ferometric observations of emission from 
ample, on the spectral energy distributions or colors or on indi- dugt around a k of 2Q ^ Q and T low . mass protostarS; 
rect measures, e.g., of the accretion rates. The most commonly afj weU as observations of a dense tracel - HCO+ 3 _ 2 , for the 
used distinction between the embedded Class and I protostars Clasf , 

I sources. By interpreting continuum data using dust radia- 

is that the former emit more than 0.5% of their luminosity at tiye tr£msfer modeK we disent a le the contrib utions from the 

wavelengths longer than 350 fim. This in turn is thought to re- envel and disks t0 the emission and estimate their masses 

fleet that they have accre ted less than half of their final mass Together whh steUar masses in f erred fr0 m Keplerian rotation in 

~ r~Z ■ T~ v t resolved line maps (see also Brinch et al., 2007; Lomme n et al.L 
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2008), we directly constrain the evolution of matter during the 
earliest embedded stages of young stellar objects. 

The formation of a circumstellar disk is clearly a key event 
in the evolution of young stellar objects, providing the route 
through which material is accreted from the large-scale enve- 
lope onto the central star. Theoretically a rotationally supported 
circumstellar disk is expected to form during the collapse of a 
dense core due to its initial angular momentum: material falling 
in from the outer regions of the core will be deflected away 
from the central star into a fiat circumstellar disk structure at the 
"centrifugal radius" where gravity i s balanced by rotation (e.g. , 
Ulrichl 119761; ICassen & Moosmanl [l98lt ITerebev etafl Il984t 
Basu. ll998l) . The centrifugal radius will grow in time, with the 
exact pace depending on the distribution of angular momentum 
in the infalling cor e, e.g., whether it can be described as solid- 
body rotation (e.g.. ITerebev et all 1 19841) or as dif ferential rota - 
tion as expected in magnetized infalling cores fe.g.. lBasa ii99a. 
Even earlier in the evolution of low-mass YSOs, an unstable 
pseudo-disk may develop as magnetic fields in the core deflect 
infalling mater ial away from the rad ial direction toward the cen- 
tral disk plane (iGalli & Shulll993allbl) . 

The existence of circumstellar disks has been used to s olve 
the so-called "luminosity problem" (e.g. lKenvon et allll990l) in 
which embedded young stellar objects are observed to be under- 
luminous compared to the expected luminosity due to release 
of gravitational energy under constant mass accretion onto the 
central star. The solution proposed by iKenvon et alj d 19901) is 
that material is accreted onto the circumstellar disk (and thus at 
larger radii with a smaller release of gravitational energy) and 
from there episodically onto the central star, e.g., related to the 
observed FU Orionis outbursts. The recent compilation of YSO 
luminosities from the Spitzer Space Telesc ope "Cores to Disks 
(c2d)" legacy program (Eva ns et~a"Q 12009) shows that even the 
most deeply embedded YSOs are under-luminous compared to 
the predictions from models - perhaps suggesting that already 
in these early phases material is not accreted directly on to the 
central protostar, but rather accumulated either in a circumstellar 
disk or accreted in another episodic manner. 

The main issue about inferring the properties of circumstellar 
disks around embedded YSOs remains to disentangle the contri- 
bution from the larger scale envelopes and circumstellar disks. 
The continuum emission from dust on the other hand provides 
an optically thin tracer, whose strength mainly depends on the 
column density and temperature of the emitting material - and 
thus potentially provides a direct probe of the dust mass of en- 
velope and/or disks if these can be separated. High angular res- 
olution line observations provide a means to complete the pic- 
ture of the structure of young stellar objects and their evolution 
by constraining the dynamical structure of the inner envelope 
and disk and through that infer the cent r al stel l ar masses (e.g. , 
Brown & Chandlert Il999l iBrinch et all 120071: iLommen et aU 



20081) . In particular, aperture synthesis observations at submil- 



limeter wavelengths are important for this goal: the higher ex- 
citation transitions of common molecules such as HCO + , HCN 
and CS observed at 1 mm and shorter wavelengt hs, probe d ensi- 
ties >10 5 cirT 3 and temperatures > 25 K (e.g.. lEvansll 19991) and 
are thus much less sensitive to the structure of the outer enve- 
lope/ambient cloud compared to ground-state and low-excitation 
transitions observable at longer wavelengths. 

For more evolved YSOs in fairly isolated regions, 
single-dish continuum observations have been used for 
comparative studies of statistica l ly significant samples of 
sources (e.gjBeckwith et allll990l:lAndre & Mont meriel,ll994t 
Andrews & Williams, 2005, 2007b). Subarcsecond resolution 



imaging using (sub)millimeter wave interfer ometers has been 
done for smaller s a mples of objects ( e.g.. | Lay et al.L [ 1994; 
i Mundv et aU 119961: IHogerheiide et all Il997bt lLoonev et all 
2000). For embedded protostars, the issue arises that both the 
disks and envelopes may contribute significantly to the emission 
within the typical single-dish beam and for such sources, inter- 
ferometric imaging is needed to address the d ust distribution on 
the siz e scales where disks may be forming. iKeene & MassonI 
(1990) demonstrated the potential of such high angular reso- 
lution millimeter observations using data from the OVRO mil- 
limeter array: by comparing the interferometric data to mod- 
els for the collapsing envelope around the L1551-IRS5 Class 
I protostar, inferred from modeling its SED, IKeene & Masson 
found that a central unresolved component was required to ex- 
plain the observed brightness profiles in the interferometric data, 
such as expected from a central circumstellar disk. In the case 
of L1551-IRS5, this compact component accounts for about 
50% of the total dust continuum emission - or if interpreted 
in the context of a circumstellar disk and envelope - a cir- 
cumstellar disk containing about 20-25% o f the total mass in 
the system. IHogerheiide et al.l dl997al 1 19981) used the OVRO 
Millimeter Array to survey the continuum and line emission at 
3 mm from a sample of 9 Class I protostellar sources. They 
showed that most of the continuum emission was associated 
with a compact, unresolved (<3") component. More recently 
these kind of studies have seen significant progress with the 
advent of n ew data and more detailed dust radiative transfer 
models (e.gjHogerheijde & SandelH 120001: lLoonev et aUl 2003 



| j0rgensen et al.1 l2004al l2005at |j0rgensenl 12004 : ISchoier et al. 
20041: lEisner et all 120051: IBrinch et all 120071: iLommen eTal 



2008). In many cases, evidence is found for similarly compact 
components - although the exact interpretation of their nature 
still depends on assumptions about the properties of the larger 
scale envelopes. Still, despite these efforts the studied samples 
have been limited to at most a handful of sources. 

In this paper, we combine continuum and line observations 
for a large sample of embedded protostars to constrain the 
masses of their three main components: the envelopes, disks and 
central stars. Specifically, we discuss single-dish and interfero- 
metric observations of a sample of 20 Class and I young stel- 
lar objects covering size scales from a few thousand down to 
a few hundred AU. We use continuum data together with de- 
tailed dust radiative transfer to constrain the masses of the larger 
scale envelopes and circumstellar disks, thereby providing con- 
straints on the evolution of the material through these embedded 
stages. In addition, we report observations of the HCO + 3-2 line 
at 267 GHz tracing dense gas in the envelopes around the Class 
I protostars - in a few cases showing evidence for Keplerian ro- 
tation in the inner envelopes/circumstellar disks. This paper is 
laid out as follows: ^introduces the selected samples of Class 
and I protostars. The details of the new SMA observations of 
the Class I sources in thi s sample are presented in $3j for the 
Class sources we refer to J0rgen sen et al.l d2007al) . $4]discusses 
the constraints on the envelope and disk properties from the con- 
tinuum data compared to generic dust radiative transfer models 
of protostellar envelope. ^5] discusses the results of the line ob- 
servations of the Class I sources - in particular, the use of HCO + 
3-2 as a tracer of Keplerian motions in the inner envelopes and 
disks around four of these sources. A discussion of outflow cav- 
ities observed in the HCO + 3-2 emission for another three of 
the Class I sources is deferred to the Appendix where also de- 
tails of the individual sources are given. Finally, §|6] compares 
the inferred masses to models for the evolution of the envelope, 
disk and stellar masses and discusses the assumptions in, e.g., 
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the assumed dust opacities, before $3 summal "i zes the main con- 
clusions of the paper. 



2. Sample and data 

2.1. Sample 

Our data consist of line and continuum observations of 10 
Class sources and 9 Class I sour ces obse r ved a s part 
of the Submillimeter Array (SMA; iHo et all |2004l3 key 
project, Protostellar Su bmillimeter Array Campaign (PROSAC; 
|j0rgensen etaT1l2007al) . 

The Class sources were observed in 3 different spectral 
setups and continuum at 1.3 mm (225 GHz) and 0.85 mm 
(345 GHz) between November 2 004 and January 2006, and the 
data for those were discussed in dJ0rgensen et al.Ll2007ah . 

The Class I YSOs were selected amongst known Class I 
sources in th e Ophiuchus and Taurus st ar fo rming regions also 
survey ed by iHogerheiid e et al.l dl997al) and Ivan Kempen et al.1 
(2009). The sample was constructed to span the expected evolu- 
tionary range o f Class I sources base d on their bolometric tem- 
peratures, 7boi dMvers&Laddl 119931) . In total eight fields were 
selected encompassing 9 YSOs previously classified as embed- 
ded Class I protostars. In addition, we include the Class I pro- 
tostar L1489-IRS, which was observed with the SMA in com- 
pact and extended configuration in the same lines as the remain- 
ing Class I source s and analyz e d usin g detailed 2D radiative 
transfer models by Brin ch etall (|2007). These surveys provide 
(sub)millimeter continuum data of 1-3" resolution with typical 
sensitivities of a few mJy beam -1 . The primary beam field of 
view of the SMA is 43-51" at 1.1-1.3 mm. 

An overview of the total sample is given in Table Q] In 
this table we give the positions of the sources from fits to 
the (sub)millimeter continuum position as well as bolometric 
temperatures and luminosities for the sources - predominantly 
based on the compilation from the Spitzer/c2d legacy program 
( Eva ns et al 1 120091) . 

In addition to the SMA data we collected additional near- 
infrared data for each Class I source: Hubble Space Telescope 
(HST) 1.6 fim images from the Near Infrared Camera and 
Multi-Object Spectrometer (NICMOS) camera were obtained 
from the HST archive at ESO for TMR1, TMC1A, IRS 43 
and GSS30. The i mages were registere d using existing near- 
infrared data (e.g., Bar sonv et al.L 11997b . We also include sin- 
gle pointing HCO + 3-2 spectra from the James Clerk Maxwell 
Telescope (JCMT ) archiv e that were originally presented by 
Hogerheii de et al.l (Il997al) and unpublished observations that 
were obtained in conne ction with the HCO + 4-3 study by 
Ivan Kempen et al.l d2009l) . For both Class and I sources, re- 
calibrated single-dish submillimeter continuum maps (Fig. [TJ 
were obtained through the JCMT/SCUBA archive legacy project 
dDi Francesco et alll2008l) . 



3. SMA observations of Class I sources 

3.1. Observations 

The observations of the Class I sources were performed between 
May 2006 and July 2007 with the Submillimeter Array (SMA) 



in two variations of its compact configuration. For all observa- 
tions, all 8 antennas were available in the array, providing pro- 
jected baselines of 9 to 64 kA (May 2006 and January 2007 ob- 
servations) and 9 to 104 kA (June and July 2007 observations). 
Typically, two sources were observed per track with the excep- 
tion of IRS 63 a nd Elias 29 that were o bserved in separate tracks 
as discussed by Lomme n et al.1 d2008l) . 

The complex gains were calibrated by observations of 1 .5- 
3 Jy quasars typically located within 15° of the targeted sources 
once every 20 minutes. The bandpass was calibrated through ob- 
servations of strong quasars and planets at the beginning and 
end of each track. The quasar fluxes were bootstrapped through 
observations of Uranus with a resulting as 20% flux calibration 
uncertainty. 

The SMA receivers and correlator was configured to ob- 
serve the lines of HCN 7 = 3^ 2 (265.886444 GHz) and 
HCO + 7 = 3^ 2 (267.557648 GHz) together with the con- 
tinuum at 1.1 mm: 1 chunk of 512 channels was centered on 
each of the lines, resulting in a spectral resolution of 0.2 MHz 
or 0.23 km s _1 . The remainder of the 2 GHz bandwidth in each 
sideband was used to record the continuum. 

The initial flagging of the data as well as bandpass, flux and 
phase/a mplitud e calibrations were performed in the Mir pack- 
agejOi, 2008) and subsequent imaging and cleaning in Miriad 
1995b . Table [3] summarizes the details of the obser- 



vations. 



3.2. Continuum data overview 



1 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics, and is funded by the Smithsonian 
Institution and the Academia Sinica. 



Fig. |2] shows the SMA continuum maps of the embedded pro- 
tostars in the sample at 1 . 1 mm for the Class I sources as well 
as the 1.3 mm con tinuum maps for the C lass sources previ- 
ously presented by J0rg ensen et al.l d2007al) . All sources, except 
IRS 54, are detected in continuum: one source, TMC1 A, shows a 
fainter secondary component, TMC1A-2 at a separation of 5.5" 
(« 770 AU) whereas the GSS 30-IRS1, -IRS2, -IRS3 system 
only shows a continuum detection for GSS30-IRS3. 

Most of the observed continuum emission has its origin in 
the inner few hundred AU of the embedded protostars. This 
is clearly demonstrated in plots of the observed visibilities as 
function of projected baseline length (Fig. [3). Most sources are 
roughly consistent with a point, or marginally resolved, source 
at the resolution of the SMA. These plots are made using cir- 
cular averages in the (u, v)-plane, i.e., implicitly assume that the 
emission is spherically symmetric. 

To estimate the structure of each source an elliptical 
Gaussian was fitted to the continuum flux in the (u, v) plane 
for baselines longer than 20 kA, where the contribution from 
the larger scale collapsing envelope becomes negligible and the 
plots of visibility amplitude vs. projected baseline length flat- 
tens (Fig. [3). The results are given in Table [3] The derived 
fluxes are slightly lower than the s ingle-dish measurements by 
lAndrews & Williams! d2005l l2007ab extrapolated to 1.1 mm as- 
suming optical thin dust with k oc \fi and y3 = 1 - but still 
more than 50% of the total flux for each source is recovered. 
This suggests that although the envelope itself contributes on 
scales smaller than the approximately 15" of the JCMT beam, 
the flux from the central compact component dominates. The 
sources that show the biggest discrepancies between the single- 
dish fluxes and interferometric fluxes (IRS 63, WL 12, GSS30 
and Elias 29) are those that show clear evidence of extended 
emission in Fig. Q] The deconvolved sizes of the continuum 
sources from the Gaussian fits suggest that they represent struc- 
tures with sizes of up to 250-300 AU diameter. 



4 J. K. J0rgensen et al.: The mass evolution of embedded protostars 

Table 1. Sample of embedded protostars. 



Source 


RA [J2000] 


DEC [J2000] 


r bo i [K] 


£*ol [Lai 


rf[pc] 






Class 








L1448-mm 


03 25 38.9 


+30 44 05.4 


77 


5.3 


220 


urr i ill. 












-IRAS2A 


03 28 55.6 


+31 14 37.1 


57 


20" 


220 


-IRAS4A-SE 


03 29 10.5 


+31 13 31.0 


43 


5.8 


220 


-IRAS4A-NW 


03 29 10.4 


+31 13 32.3 


- 


- 


220 


-IRAS4B 


03 29 12.0 


+31 13 08.1 


55 


3.8* 


220 


-IRAS4B' 


03 29 12.8 


+31 13 07.0 


— 


— 


220 


LI 527 


04 39 53.9 


+26 03 09.8 


59 


2.0 


140 


L483 


18 17 29.9 


-04 39 39.6 


60 


9.0 


200 


B335 


19 37 00.9 


+07 34 09.7 


60 


3.0 


250 


L1157 


20 39 06.3 


+68 02 16.0 


42 


6.0 


325 






Class 1 








L1489-IRS 


04 04 42.9 


+26 18 56.3 


238 


3.7 


140 


TMR1 


04 39 13.9 


+25 53 20.6 


144 


3.7 


140 


TMC1A 


04 39 35.2 


+25 41 44.4 


172 


2.2 


140 


GSS30-IRS3 


16 26 21.7 


-24 22 50.6 


300 


14'' 


125 


GSS30-IRS1 


16 26 22.2 


-24 23 01.9 






125 


WL 12 


16 26 44.2 


-24 34 48.7 


440 


2.6 


125 


Elias 29 


16 27 09.4 


-24 37 20.0 


350 


41 


125 


IRS 43 


16 27 26.9 


-24 40 50.6 


310 


6.0 


125 


IRS 54 


16 27 51.8 


-24 31 45.4 


740 


2.5 


125 


IRS 63 


16 31 35.7 


-24 01 29.6 


530 


3.3 


125 



"For NGC 1333-IRAS2A lEvans et alj d2009h quote a luminosity of 76 L Q . Estimating the luminosity from scaling the 70 /im flux using the 
approach of Dunham et al. (2008D gives a luminosity of 19.6 L Q , which agrees well with previous estimates from modeling the SED in the range 
16-20 L Q . The higher luminosity in the c2d compilation is a result of the SED of IRAS2A being unconstrained at 24 /im (M. Dunham, priv. 
comm) - and we therefore adopt a luminosity of 20 L Q for this source. *The Spitzer maps of NGC 1333-IRAS4B shows significant emission 
from a nearby outflow shock within the typical single-dish beam. As total luminosity for th is source, we adopt the sum of the two components, 
SSTc2d J032912.06+311305.4 and SSTc2d J032912.06+311301.7, from lEvans etal] J20091) . The emission from GSS30-IRS3 (LFAM1) and 
GSS30-IRS1 are not separated in the Spitzer catalog (the latter dominates the emission and its luminosity is quoted here). 

Table 2. Log of observations of Class I sources. 



Sources 



Date 



Gain calibrator (flux) 



Beam size Continuum RMS 



T" 



IRS 63 2006 May 15 J1626-298 (1.3 Jy), J1517-243 (2.4 Jy) 4.0"x2.3" 5.5 mJy beam 

Elias 29 2006 May 17 J1626-298 (1.7 Jy), J1517-243 (3.2 Jy) 4.0"x2.3" 3.6 mJy beam 

TMR1, TMC1A 2007Jan02 J0530+135 (2.8 Jy), 3clll (2.3 Jy) 2.7"x2.4" 2.6 mJy beam 

IRS 43, IRS 54 2007 Jun 27 J1626-298 (1.5 Jy), J1517-243 (1.6 Jy) 3.9"x2.0" 2.8 mJy beam 

WL 12, GSS30 2007 Jul 04 J1626-298 (1.3 Jy), J1517-243 (1.4 Jy) 2.8"x2.7" 2.6 mJy beam 



3.3. Line data overview 

Fig. [4] shows an overview of the HCO + 3-2 spectra toward 
the continuum positions of each source from the interfero- 
metric observations, as well as single-dish spectra from the 
JCMT archive. The HCO + 3-2 line profiles vary significantly 
in strength between different sources, but also relatively be- 
tween the single-dish and interferometric data. In the single-dish 
data all sources show HCO + 3-2 lines with intensities of 1.5- 
1 K km s~ 1 , but in the interferometer data only 6 of the 9 Class I 
sources have clearly detected lines. IRS 43 shows the highest 
intensity and broadest line profile in the HCO + SMA obser- 
vations. Interestingly, the interferometric spectra toward TMR1 
and TMC1 A only probe the blue-shifted emission relative to the 
single-dish systemic velocity. Fig. [5] shows SMA maps of the 
HCO + and HCN emission integrated over 4 velocity intervals 
around the systemic velocity. Again, clear differences are found 
between the sources. IRS 63 and IRS 43 show red- and blue- 
shifted emission around the continuum position whereas TMR1, 



TMC1A and GSS30-IRS1 display emission extended from one 
side relative to the source position. Toward Elias 29 a com- 
plex structure with large-scale extended HCO + emission from 
the nearby ridge is seen. For that source material with the most 
extreme red/blue-shifted velocities around the continuum posi- 
tion trace the Keple rian rotation in the inner envelope/disk (see 
iLommenet all 120081) . 



4. Analysis of continuum data: disk vs. envelope 
masses 

With large continuum surveys at hand the first step is to dis- 
entangle the envelope and disk contributions to the dust con- 
tinuum emission and from this derive the disk and envelope 
masses. Ideally fully self-consistent models can be derived for 
all sources taking into account the full spectra l energy distri- 
bution s as well the in t erfero metric data (e.g., J0rgensen et all 
2005a; Lomm en et aU [2008). However, much information can 
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Fig. 1. Spitzer Space Telescope 3.6, 4.5 and 8.0 fim images (blue, green and red) of each of our sources with SCUBA 850 yum emis- 
sion ( |Pi Francesco et aijj2008l) shown as contours. The contours are given in steps of 0. 1 Jy beam 1 - except for L1448, IRAS2A, 
and IRAS4A/B where the (d otted) contours a r e given in step s of 0.3 Jy beam" 1 . For furth er details about Spitzer observ ations of indi- 
vidual sources see: Pers eus: |j0rgensen et al] d2006l l2007b|) andlGutermuth et all (120081). Taurus : iPadgett et al] (120081), Ophiuchus: 
lEvans et al.l (120091) and Ifergensen et al] d2008l) . L1527: iTobin et al] (120081) . B335: IStutz etaTI (|20()8|), and L1157: lLoonev et al] 



be derived from just simple comparisons of the single-dish and 
interferometric fluxes. 

As discussed by, e.g., iTerebey et al j (1 19931) the continuum 
emission from a cloud core or circumstellar disk is well approx- 
imated by low optical depth in which case the intensity pro- 
file is simply found by integrating the temperature and density 
pr ofile of the core or disk along the line of sight (e.g., Eq. 1 
of ITerebey et al.l fl993l) . The complication naturally arises that 



both the density and temperature vary throughout protostellar 
disks and envelopes. Dust radiative transfer models calculate the 
dust temperature distribution self-consistently and furthermore 
produce images at different wavelengths, which can be directly 
compared to the observations. 

In this section we use such dust radiative transfer models to 
derive the masses of the envelopes and disks around the sam- 
ple of Class and I protostars. We first explore one-dimensional 
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10 5 -5-10 10 5 -5-10 10 5 -5-10 10 5 -5-10 10 5 -5-10 10 5 -5-10 

RA offset ["] RA offset ["] 



Fig. 2. Overview of the 1 .3 mm continuum emission from the Class sources (left group of panels; see also J0rgensen et al. 2 007af) 
and 1.1 mm continuum emission from the sample of Class I sources (right group of panels; this paper). Contours are given from 3<x 
to 18cr in steps of 3<x (black contours) and onwards from there in steps of 6<x (white contours). In the panels for IRS 54 and GSS30 
the Spitzer positions for the known YSOs are indicated with the plus-signs. 

Table 3. Results of elliptical Gaussian and point source fits to the continuum visibilities for the Class I sources. 



Source Position Size", FWHM Flux (density) 

RA DEC ^l.lmra F 1.1 mm S 850 /an(15")' 

(J2000.0) (J2000.0) H [mJy] [mJybeanr 1 ] 



TMR1 


04 39 13.91 


+25 53 20.63 


d 


79" 


480 


TMC1A 


04 39 35.20 


+25 41 44.35 


0.69x0.40 


256 (246) 


780 


TMC1A-2 


04 39 34.80 


+25 41 43.87 


3.7x1.1 


48 (26) 




IRS 43 


16 27 26.91 


-24 40 50.62 


1.9x0.13 


75 (56) 


990 


IRS 54 


16 27 51. 80* 


-24 31 45. 4 d 








IRS 63 


16 31 35.65 


-24 01 29.55 


0.63x0.0096 


473 (459) 


1220 


GSS30-IRS3 


16 26 21.71 


-24 22 50.63 


1.4x1.3 


204 (144) 


990 


WL 12 


16 26 44.19 


-24 34 48.74 


1.9x1.5 


104 (64) 


360 


Elias 29 


16 27 09.44 


-24 37 19.99 


5.8x1.8 


109 (48) 


590 



"Deconvolved size from Gaussian fit. *Flux density from Gaussian fit to SMA data in the (u, v)-plane. Number in parentheses from point source fit. 
r Peak flux at 850 yum in ICMT 15" single-dish beam. f 'TMRl unresolved in SMA observations. Only result from point source fit listed. IRS 54 
undetected at the SMA; position from Spitzer Space Telescope observations. 



dust radiative transfer models of protostellar envelopes: we com- 
pare the emission from a spherically symmetric envelope around 
a central source of heating (i.e., protostar) at the scales observed 
with single-dish telescopes compared to our interferometric ob- 
servations ( 94.lt and furthermore use the models to establish 
a direct relationship between the single-dish submillimeter flux 
and mass of an envelope around a given protostar as function of 
its luminosity and distance ( §4.2t . Finally, we use these results 
to separate the envelope and disk emission from the envelopes 
and disks and estimate their masses for each of the objects in 
our sample ( fl4~3l and jp~4l ). 



4.1. Dust radiative transfer models; how much does an 
envelope contribute to the compact emission? 

First, to test the hypothesis that most of the flux in the inter- 
ferometric data at baselines of 50 kA is in fact from the central 
compact component we calculated a set of self-consistent mod- 
els for the dust in protostellar env elopes using the Dusty radiative 
transfer code dlvezic et al.L[l99 9). Dusty solves the temperature 
profile for a spherically symmetric envelope around a central 
source of heating given the power-law envelope density profile, 
input spectrum of the central heating source and dust proper- 
ties. The code furthermore calculates images at user-specified 
wavelengths, which can be compared to observed images to iter- 
atively place constraints on, e.g., the density profile. It has pre- 
viously been used to constrain the physical structures of the en- 
velopes around a large sample of deeply embedded protostars 
(j0i-g ensenetan.l2OO2l) . 
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Fig. 4. HCO + 3-2 single-dish spectra from the JCMT archive (left) and toward the continuum positions of each source from the 
SMA observations (right). A Gaussian was fitted to each single-dish spectrum for each source and the centroid velocity as well as 
velocities +2 km s _1 from this are shown with dotted lines in each panel and used to calculate moment maps in Fig. [5] 
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Fig. 3. Plot of the continuum flux as function of projected base- 
line length. The dotted histograms indicate the zero-amplitude 
signal, i.e., the anticipated signal in the absence of source emis- 
sion. The open circles indicate the single-dish peak flux [in 
Jy beam 1 ] from the JCMT/SCUBA maps toward each source 
extrapolated to 1 . 1 mm. 



We fixed the effective temperature of the central black-body 
to 1500 K mimicking the combined contribution from a low ef- 
fective temperature young star and a disk. Since most of this 
emission is reprocessed in the larger-scale cold envelope, the ex- 
act shape and temperature of the input spectrum does not af- 
fect the distribution and strength of the emerging submillimeter 
emission (see also lJ0rgensen et al.1 120021 l2005bb ISchoier et al.l 
l2002t IShirley etaf]l2002h . We furthermore fixed the outer en- 
velope radius to 8000 AU (typical of envelopes in clustered 
regions such as Ophiuchus and Perseus) and adopted enve- 
lope density profiles with n cc r~ p with p = 1.5. Such den- 
sity profiles are expected for free-falling envelopes, e.g., within 
th e col l apse e xpansion wave in the inside-out collapse model 
of IShul (1 19771) . Detailed models of envelopes of Class and I 
sources constrained by SCUBA mappin g observations find den- 
sity profiles with p = 1.5 - 2.0 (e.g., iJflrgensen et all [2002: 
IShirlev et all 120021; lYoung et all 12003b . Those models do not 
include the contributions from a central disk, which can steepen 
the derived density profile index artificially by 0.2-0.5. A power- 
law index of 1.5-1.8 for the envelope density profile therefore 
appears to be consistent with the envelope structure at these 
scales. As in previous work, we adopt the dust opacity law of 
lOssen kopf & Hennin g| (11994 for grains with thin ice mantles, 
coagulated at a density of «h, ~ 10 6 cm~ 3 . This dust opacity 
law has k^q^ = 0.0182 cm 2 g -1 (dust+gas) and scales approx- 
imately with frequency as k v oc v lJ at 1 mm. Since we for now 
just deal with the fluxes at submillimeter wavelengths, the de- 
rived masses can be directly rescaled using a different dust opac- 
ity law. 

With these fixed parameters we then calculated a set of dust 
radiative transfer models varying the inner radius, luminosity 
and mass of the envelope. We tested models with inner radii, 
Ri = 10, 25, 50, 100, 200 AU, envelope masses between 0.01 and 
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Fig. 5. Overview of the HCO + 3-2 (left) and HCN 3-2 ( right) emission from the sample of Class I sources. The contours (shown in 
3<t intervals) indicate emission integrated over intervals from -4 to -2 km s _1 (blue), -2 to km s _1 (green), to +2 km s _1 (orange) 
and +2 to +4 km s -1 (red) relative to the systemic velocities from the Gaussian fits to the single-dish spectra. 



0.5 M Q and luminosities of 1 and 5 L©. Together these parame- 
ter ranges cover models representative for the majority of our 
sources. 

For each model, the radiative transfer code calculates the 
temperature profile self-consistently as well as images at the 
wavelengths of our single-dish and interferometric observations 
(850 /im and 1.1 mm, respectively). The images at the wave- 
length of our interferometric data were subsequently multiplied 
by the SMA primary beam and Fourier transformed, mimick- 
ing the interferometric observations. From these observations 
we extracted the flux at baselines of 50 kA (3.5" diameter) as 
an estimate of the contribution to the SMA 1 . 1 mm flux from 
the envelope. Likewise, images at 850 yum were convolved with 
the SCUBA beam and the peak flux in the central beam was ex- 
tracted as an estimate of its single-dish 850 /-im peak flux. 

Fig. [6] shows the ratio of the 1.1 mm interferometric fluxes 
at baselines of 50 kA and single-dish peak fluxes at 850 fim as 
function of single-dish peak flux. As shown, the contribution of 
the envelope to 1 . 1 mm emission at baselines of 50 kA is at most 
only 4% of the modeled 850 //m peak flux. If the envelope is 
steeper, n oc r~ L8 , which seems to be the steepest envelope struc- 
ture supported by the modeling of the single-dish data, the enve- 
lope contribution to the 1 . 1 mm emission at baselines of 50 kA 
increases to 8% of the peak single-dish flux at 850 fim. 

As expected models with large inner radii resolved by the 
interferometer show less compact flux lowering the ratio of the 
interferometer-to-single-dish flux. This increase in flux with de- 
creasing inner radius reaches its maximum of 4% for an inner ra- 
dius of about 25 AU. Material at smaller inner radii is not probed 
by the existing interferometric observations: this reflects the dis- 
tribution of mass in the power-law envelope ensuring that most 
mass is at large scales, or that most of the dust continuum emis- 
sion is due to mass located at the typical size scale probed by the 
interferometer. 



Using these results we can estimate an upper limit to the en- 
velope contribution at 50 kA for the observed sources by taking 
4% of the 850 fim JCMT peak flux and comparing that directly 
to the observed 1.1 mm interferometric flux at 50 kA: the ra- 
tio between these two numbers then gives the relative envelope 
contribution to the 1.1 mm interferometric data. This compari- 
son confirms that the 1 . 1 mm SMA measurements of the Class I 
sources (Table [3} predominantly contain emission from the cen- 
tral disk: in median, 87% of the observed interferometer flux is 
from the disk, as also suggested by their flat visibility amplitudes 
as function of projected baseline length in Fig. [3] For compari- 
son, the Class sources have a median disk contribution of 68% 
of their observed interferometer flux. The results in Fig. [6] are 
shown for a distance of 125 pc. The ratio between the interfer- 
ometer and single-dish peak flux measures the relative surface 
brightnesses on the spatial scales probed by each and does not 
change significantly when the distance increases - as long as the 
inner envelope cavity is not large enough to be resolved (Fig.[7]l. 



4.2. Calibration of envelope mass - continuum flux relation 

With a separation of the envelope and disk contribution to the 
single-dish and interferometric observations, the next task is to 
derive masses from each component using the observed fluxes. 
For prestellar cores it is usually a good approximation to assume 
that their dust is optically thin and isothermal with a temper- 
ature of 10-15 K and from that estimate the mass of the core 
from single-dish observations of their submillimeter continuum 
fluxes. However, an envelope with an internal heating source, 
such as a central protostar, will have a distribution of tempera- 
tures - which will strongly depend on the source luminosity. 

A way to circumvent these problems is again to use the re- 
sults from the detailed dust radiative transfer models: Fig.[8]com- 
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Fig. 6. Ratio of interferometric fluxes (1.1 mm flux at 50 kA; 
S^ofcU.i mm) and single-dish fluxes (850 fim peak flux in 15" 
SCUBA beam; S 15 » 850A(m ) as function of single-dish flux for a 
set of models varying the inner radius (10 AU, 25 AU, 50 AU, 
100 AU and 200 AU shown with red, orange, green, light blue 
and dark blue lines, respectively; see also upper panel), source 
luminosity (1 and 5 L Q shown with thin and thick lines, respec- 
tively) and envelope masses (along each model line). Envelopes 
have density profiles, n oc r~ p with p = 1 .5 ( upper) and p — 1 .8 
(lower) and masses ranging from 0.01 M at low S is« to 0.5 M 
at high S i5". 



pares the relations between the envelope mass, peak flux and lu- 
minosity for the protostellar envelopes determined from the ra- 
diative transfer modeling. As shown these dust radiative transfer 
models give rise to simple power-law relationships between the 
envelope mass and observables which can be expressed as: 



0.44M o 



-0.36 



15 ,850/jm 



1 Jy beam 



-l 



125 pc 



1.2 



(1) 



Together with the luminosity, the above expression makes it pos- 
sible to evaluate the mass of each protostellar envelope from its 
peak flux only - once the contribution from the disk emission to 
the single-dish flux has been subtracted (see also Sect. 14. 3b . This 
relation is particularly useful in crowded regions where it is pos- 
sible to estimate the peak flux of a given source, but where it is 
difficult to evaluate its outer radius and thus its extended flux due 
to confusion. The luminosity can either be derived from the full 



Fig. 7. Ratio between the interferometer flux at 1.1 mm and 
single-dish beam peak flux at 850 pm as function of distance 
to the source for an envelope with p = 1.5, R oM = 8000 AU, 
R, = 25 AU (dashed line) and 100 AU (solid line) and an enve- 
lope mass w 0.3 M Q . 



spectral energy distribution of each source or by using the em- 
pirical correlat ion between source in ternal luminosity and 70 pm 
Spitzer flux bv lDunham et al.l (120081) . 

These results are generally in good agreement with the mod- 
els discussed by iTerebev et al.l d 1 993b : they find a relationship 
between the flux and luminosity and distance which for a given 
envelope mass scales roughly as F oc L°- 2 d^ 0J compared to 
our results which for a given envelope mass are equivalent to 
F oc L 03 d~ h0 , The different slope with respect to distance is 
caused by the sligh tly flatter density profile predicted by the 
ITerebev et alj (1 19841) model: for a power-law envelope similar 



to ours, Terebev et al.l (Q993) find that F oc d~ 09 . The remain- 
ing difference in the exponent of this relation as well as the 
difference in the exponents for the luminosity is caused by the 
envelopes being slightly optically thick to their own radiation 
which increases the temperature (and thus emission) on small 
scales and causes the underlying temperature profile to depart 
from the single power-law applicable for an envelope fully opti- 
cally thin to its own radiation. 

The derived expressions are particularly useful for interpret- 
ing large samples of sources detected, e.g., through the large- 
scale submillimeter maps provided by ongoing and upcoming 
surveys with bolometer arrays such as LABOCA and SCUBA 2. 
One of the main goals of those surveys are naturally to deter- 
mine core masses and this is often done using single isother- 
mal cores with two temperatures: one for cores with and one 
without embedded proto stellar sources (e.g.. lEnochetalll2007t 
lHatchell & Fulieil I2008T) . However, with Eq.Q] we can take the 
dependence of the source submillimeter flux on the luminosity 
into account: according to Eq.Q]a factor 2.3 decrease in mass is 
expected due to the increase of the luminosity by a factor 10 for 
the same peak flux - or put in another way: if the mass for a 1 L Q 
core is derived using a temperature of 15 K, that for a 10 L G pro- 
tostellar core should adopt a temperature of 23.5 K. Assuming 
a similar temperature for all protostellar cores would therefore 
lead to a mass-luminosity relation, L ~ M 2 8 . This is slightly 
steeper but close to the relation found by Hatch ell et al.l ([2007 ) 
with L ~ M 1 96 with an uncertainty of +0.36 in the power-law 
index of the observationally derived relation. 
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Fig. 9. Single-dish and interferometer fluxes as well as disk and envelope masses as function of source bolometric temperatures. 
The upper panels show the ratios of the interferometer and single-dish flux (left) and the interferometer flux (right) and the lower 
panels the ratio of the disk and envelope mass (left) and disk mass (right). Class sources are shown with red symbols and Class I 
with blue symbols. In the lower left panel, L483 (Class source with envelope-only emission; see discussion in text) is shown with 
the red triangle. In the lower right panel, IRS 54 (Class I source with no detected interferometric or single dish emission) is shown 
with a 3 cr upper limit to its disk mass (0.001 M Q ) from the SMA observations with the blue triangle. 



4.3. Determining envelope and disk masses from continuum 
measurements 

With the resuts from the previous section we can now compare 
the derived compact and single-dish fluxes across the sample. 
The Class and I sources were not observed at exactly the 
same wavelengths at the SMA. The Class sources were ob- 
served at 1.3 mm (230 GHz) and 0.8 mm (345 GHz) whereas 
the Class I sources discussed here were observed at 1.1 mm 
(270 GHz). Class and I sources have been observed at 850 fim 
with SCUBA on the JCMT. To take this difference into account 
we interpolate the Class poi nt source fluxes estimate d on base- 
lines longer than 40 k/l from J0r gensen et al.l d2007a) using the 
derived spectral slope of 2.5 from that paper, i.e., less steep than 
for the envelope. With this correction the fluxes of each source 
at similar angular scales are compared directly (Fig.|9]l. 

To estimate the masses of both components, a few more 
steps are required: both the envelope and disk contribute to the 
observed emission in both the single-dish and interferometer 
beams. However, with the above results we can quantify the 
maximum contribution from the envelope on the 1 . 1 mm inter- 
ferometric flux at 50 k/l as no more than about 4% of the enve- 



lope single-dish peak flux at 850 /jm. Using this upper limit, we 
can write the two equations with the envelope flux at 850 /im, 
Serw, an d disk flux at 1.1 mm, 5 disk as two unknowns that can 
easily be solved for: 

S 50U = S disk + C ■ S env (2) 
V = Sd,sk (1.1/0.85)° + S ellv (3) 

In these equations a is the spectral index0 of the dust contin- 
uum emission at millimeter wavelengths in the disks found to be 
» 2.5 from the observati ons of the Class sources in our sample 
( J0rg ensen et al.Ll2007al) . c is the fraction of the envelope 850 yum 
single-dish peak flux observed at 1 . 1 mm by the interferometer 
at 50 k/l. For this discussion we adopt the upper limit to c = 0.04 
from the above discussion for p = 1.5. We note that this contri- 
bution generally is small enough that the direct measurement by 
the interferometer is a good first estimate of the actual compact 
flux not associated with the infalling (n oc r~ L5 ) envelope. 

To derive envelope masses, a few additional effects need to 
be taken into account here: as shown above the envelope masses 

2 Not to be confused with the infrared slope, am, measured from the 
Spitzer data at 2.2 to 24 pm. 
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Fig. 8. Results of dust radiative transfer calculations for simple 
envelope models at d — 125 pc. Upper: the relation between en- 
velope mass and peak flux in the JCMT beam for envelope mod- 
els with 0.5, 1, 2, 3 and 5 L Q from left (red) to right (blue): each 
point corresponds to a model in our grid with increasing line-of- 
sight optical depths at 100 pm ranging from 0.01 to 0.20 in steps 
of 0.01, from 0.20 to 0.50 in steps of 0.05 and at 1.0. Overplotted 
also are power-law fits with M env oc S L2 „ . Lower: Relation 

r 15 ,850/jm 

between the M em /S 1 2 „ normalization factor and luminos- 

15 ,850/jra 

ity. Overplotted is the best fit M env /S ^ q oc L®^ 6 power-law 
relation between these quantities. 



scale with distance as d xl for a fixed peak flux, whereas the 
disk masses (largely unresolved) are expected to scale as d 2 . 
In addition, since the submillimeter continuum emission mea- 
sures the dust mass weighted by temperature, the derived enve- 
lope masses are found to depend on the luminosity (see, Eq. []]). 
Taking this into account, however, it is possible to directly esti- 
mate the masses of the envelopes from the observed fluxes once 
the contribution from the disk has been estimated through Eqs.|2] 
and[3](Table|4|l. For the disk masses we assume that the dust con- 
tinuum emission is optically thin and coming from an isothermal 
disk with a temperature of 30 K and a dust opacity at 1 . 1 mm 
similar to that of the dust in the protostellar envelopes. We re- 
turn a discussion of these assumptions in §16.11 



Table 4. Envelope and disk masses for the full sample of Class 
and I sources. 



Source 


M " 


Mdisk* 


M staI c 




[M ] 


[M Q ] 


[M Q ] 




Clas, 







L1448-mm 


0.96 


0.044 




1 [\AOZn 


1 n 
1 .u 


U.UJU 




IRAS4A1 


4.5 d 


0.46 




IRAS4A2 


_d 


0.25 




IRAS4B 


2.9 d 


0.24 




IRAS4B' 


_d 


0.089 




L1527 


0.56 


0.029 




L483 


0.82 


0.0 e 




B335 


0.93 


0.017 




LI 157 


0.94 


0.12 






Class I 




T 1 A OH TD C 


0.11 


0.018 


1 ■if 


TMR1 


0.11 


0.0094 




TMC1A 


0.12 


0.035 


0.35-0.7* 


GSS30-IRS3 


0.096 


0.026 




WL 12 


0.054 


0.011 




Elias 29 


0.047 


0.011 


2.5* 


IRS 43 


0.026 


0.0081 


1.0' 


IRS 54^ 




<0.001 




IRS 63 


0.098 


0.053 


0.36'' 



"Envelope mass determined from the direct scaling of the ICMT peak 
flux (minus the disk contribution). *Disk mass assuming a uniform 
dust temperature of 30 K. ^Central mass from Keplerian velocity pro- 
files (see specific references for details). rf IRAS4Al and IRAS4A2, re- 
spectively IRAS4B and IRAS4B' have larger scale circumbinary en- 
velopes. The compact flux in L483 is consiste nt with being solel y 
from the larger scale envelope. ^Stellar mass from B rinch et all 02007). 
envelope and disk masses calculated under same assumptions as 
for the rem a inder of the sources usi ng the continuum data from 
iBrinch et all ferown & Chandled dl999h . iLommenetaD d2008h . 'This 
work. TRS 54 is undetected in both interferometric and single-dish 
maps. A 3cr upper limit to its disk mass is estimated from the SMA 
observations. 



4.4. Resulting disk and envelope masses 

Fig.|9]compares the interferometric and single-dish fluxes scaled 
to 1 . 1 mm as well as disk and envelope masses for all the sources 
in our sample. As expected the deeply embedded, Class 0, 
protostars have larger fluxes in the single-dish relative to the 
smaller interferometric beam, reflected in a lower Ssou/SiS" 
rat io (upper left panel of Fig. |5J- A similar effect was noted 
by [Ward-Thompsonf (|2007l) who compared interferometric and 
single-dish fluxes at 850 pm from literature studies of a sample 
of 9 YSOs (6 Class 0, 2 Class I and 1 Class II sources). Taken by 
itself this can naturally reflect one or two effects: ( i) the envelope 
masses for these sources are larger also resulting in their redder 
SEDs or ( ii) the disk masses are smaller for the more deeply em- 
bedded protostars (e.g., if significant disk growth occurs from 
the Class through I stages). The interferometric fluxes do not 
show such an increase going from the Class to I stages, how- 
ever (upper right panel of Fig. [5). This argues against a signifi- 
cant increase in disk masses being the explanation. 

The lower panels of Fig. [9] compare the envelope and disk 
masses for all the sources in the Class and I samples as a func- 
tion of bolometric temperature. Again, the disk/envelope mass 
ratio is seen to be lower for the more deeply embedded Class 
protostars reflecting their higher envelope masses at comparable 
disk masses as those in the Class I stages (lower left panel of 
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Fig. Typically the Class protostars have disks that contain 
about 1-10% of the total disk+envelope mass whereas the disks 
around the Class I sources contain 20-60% of the disk+envelope 
mass. In both stages most of the derived disk masses are in the 
range of ~0.01 to 0.1 M Q with the median Class disk mass of 
0.089 M and the median Class I disk mass of 0.011 M Q . The 
drop in envelope masses from the Class to I stages and the rel- 
atively low disk masses, suggests that a significant fraction of 
the mass of the central star is attained already during the most 
deeply embedded stages. 

The results above do not show evidence for an increasing 
disk mass going from the Class sources to the Class I sources 
in the observed samples. Thus, either the disks build up rapidly 
during the early, deeply embedded stages or their physical prop- 
erties (e.g., grain properties, temperatures ) are significantly dif- 
ferent than tho se in more evolved stages. lAndrews & Williams! 
(2005, 2007b) examined the relation between the disk masses 
for a large sample of Class I and II sources in Ophiuchus and 
Taurus and found evidence for a decrease of the disk masses by 
a factor 3 between these evolutionary stages from a median disk 
mass of 0.015 M Q to 0.005 M . 

Although most sources show compact emission, there are ex- 
ceptions: the interferometric flux of the L483 Class source is 
found to be consistent with envelope-only emission and does not 
require a central compac t source. This is similar to the conclu- 
sion of J0rgensen (2004) who modeled 3 mm continuum emis- 
sion of the source from OVRO observations. It is also in agree- 
ment with the measured spectr al slope of the SMA c ontinuum 
flux between 1.3 and 0.8 mm ( J0rge nsen et al.L l2007al) . which 
shows an index of * 4, consistent with optically thin dust in a 
larger scale envelope, contrasting the other Class sources in 
the sample which have lower spect ral indices in the range of 2- 
3 (Fig. 3 of |j0rgensen et al.ll2007al) . It is puzzling that no direct 
evidence of a dust disk is seen for this source, given that it has 
a well-established protostellar outflow - which even seems to 
be in the process of actively dispersing the envelope. This may 
be a case, however, where we overestimate the contribution by 
the envelope on the longest baselines from the radiative trans- 
fer models. If the envelope, for example, has a larger flattened 
inner region or cavity, the envelope contribution on the longer 
baselines would be smaller and a central circumstellar disk re- 
quired to explain the submillimeter emission. The mass of this 
disk would in that case be <0.01 M Q - derived from the long 
baseline interferometric data with out taking the larger sc ale en- 
velope contribution into account dj0rgensen et al.Ll2007al) . 

In a recent study iGirart et all d2009i) modeled SMA contin- 
uum data for the Class protostellar binary L723 and likewise 
found that a disk was not required to explain the compact con- 
tinuum emission. Scaling their observed flux at 1.35 mm from 
the SMA at baselines of about 50 kA to our reference wave- 
length at 1.1 mm gives a flux of about 35 mJy - or about 3- 
4% of the observed SCUBA peak flux of 1 . 1 Jy beam " 1 from 
the JCMT/SCUBA legacy catalog dDi Francesco et all 120081) . 
Within our framework this fraction is exactly what is expected 
for envelope only contribution at these baselines , and in agree- 
ment with the interpretation bv lGirart et alj d2009b that no disk is 
required to explain the SMA continuum observations, but with 
a note that if the envelope would have an inner flattened region, 
an additional central dust component would be required - again 
with a typical mass < a fewxO.01 M Q . 



5. Analysis of line data: stellar masses 

In the previous section masses were determined for two of the 
three YSO components - their disks and envelopes. The resolved 
line observations obtained for the Class I sources provide infor- 
mation about the dynamical structure of their inner envelopes 
and disks and thereby a handle on the mass of the third compo- 
nent, the central star. 

As mentioned above, significant differences exist between 
the components traced by the line observations for the different 
sources in the sample. Fig.QXjcompares the near-infrared 1.6 /mi 
images to the HCO + emission in TMR1, TMC1A, GSS30-IRS1 
and IRS 43. The first three of these sources show a clear align- 
ment with the near-infrared scattering nebulosities whereas the 
fourth, IRS 43, shows the HCO + (and also HCN) emission 
aligned with the dark lane across the source. Likewise, in the first 
order moment (velocity) maps (Fig. ITTb. TMR1, TMC1A and 
GSS30 show significantly different morphologies from IRS 43 
and IRS 63: the former sources show evidence of blue-shifted 
emission over most of their maps, with velocities decreasing to- 
ward the systemic velocities at the largest distances. In contrast 
IRS 43, in particular, shows a clear large-scale velocity gradi- 
ent relative to the continuum peak position and systemic veloc- 
ity, potent ially indicative of rotation as also inferr ed for IRS 63 , 
Elias 29 dLommen et all 120081) and L1489-IRS dBrinch et all 
12007b . In the following we focus on the use of HCO + as a tracer 
of rotation in these sources and include a small discussion of the 
sources showing outflow cavities and the limitations of HCO + as 
a tracer in AppendixlAl Some details about the individual Class I 
sources are given in AppendixlBl 

5.1. Keplehan rotation and stellar mass in IRS 43 

Four sources in our sample, IRS 43, IRS 63 a nd Elias 29 
(Fig.© as well as L1489-IRS dBrinch et all I2007I) . show elon- 
gated HCO + emis sion stretching acro ss t he location of th e YSO . 
As discussed by iBrinch et all d2007l) and iLommen et al! I (I2008I) 
the velocity gradients toward L1489-IRS, IRS 63 and Elias 29 
can be interpreted as having their origin in the inner enve- 
lope/circumstellar disk around each source. A number of effects 
point toward this also being the case for IRS 43: 

1. Both HCN and HCO + show an elongated structure with 
a large scale velocity gradient around the continuum 
peak/systemic velocity. 

2. A Gaussian fit to the continuum data in the (m, v) plane shows 
a structure which is elongated in the same direction with a 
position angle of -70° (measured from north toward east): its 
deconvolved major axis is about 2" (280 AU) and its minor 
axis about a tenth of an arcsecond, consistent with emission 
from a disk seen at a high inclination angle, clo se to edge-on . 

3. The direction of the Herbig-Haro objects dGrosso et all 
1200 lb and proposed radio thermal jet (Gira rt et al.Ll2000h are 
with a position angle of 20-25°, perpendicular to this ex- 
tended structure (Figs.[TT1and |B.l| ). 

The narrow continuum structure and the alignment with 
the HCO + 3-2 emi ssion is also seen in the IRS 63 data from 
ILommen et al.l d2008l) . For that source, the velocity field could 
be well-fit with a Keplerian profile around a 0.37 M Q central ob- 
ject for an inclination of 30° (as hinted by modeling of the spec- 
tral energy distribution of the source). A position-velocity plot of 
the HCO + 3-2 emission toward IRS 43 is shown in Fig. [12] The 
position velocity diagram shown here has been extracted along 
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Fig. 11. Moment one (velocity) maps of the HCO + emission for IRS 43, IRS 63, TMR1, TMC1A and GSS30-IRS1. In all panels 
the plus-signs indicate the location of the continuum positions (for GSS30-IRS1 the near-/mid-infrared source). In the IRS 43 and 
IRS 63 panels the white line shows the direction of the continu um structure. In th e IRS 43 panel, t he black arrow furt hermore shows 
the direction of the embedded near-IR Herbig-Haro objects dGrosso et all 1200 ll) and thermal jet dGirart et all uOOO). 



the major axis of the HCO + emission at a position angle of -70°. 
The velocity profile is consistent with Keplerian rotation. 

The extent of the HCO + emission (« 15") suggests that its 
origin is not purely in the central disk, but also the inner regions 
of the rotating envelope (r < 1000 AU). This is further supported 
by the comparison between the single-dish spectrum and a spec- 
trum extracted from the interferometric data integrated over the 
extent of the JCMT single-dish beam (Fig. [13). About 40% of 
the single-dish flux is recovered by the interferometer, suggest- 
ing the presence of some extended emission resolved out by the 
interferometer - although this fraction is significantly less than 
typically s een, e.g., in the more de eply embedded Class proto- 
stars (e.g., J0rg ensen et alll2007al) . 

To derive the mass of the central object a position-velocity 
curve is extracted by fitting the position of the peak emission for 
each channel in the (m, v)-plane at baselines larger than 20 kA 
and projecting this on the major axis of the HCO + emission 
(right panel, Fig.[T2l. To this curve, a^- 2 -fit is performed with a 
straightforward Keplerian rotation curve with two free parame- 
ters, the systemic velocity and central mass. A best fit is obtained 
for a central mass of 1.0 M Q (2<x confidence levels of +0.2 M ) 
and a systemic velocity of 4.1 km s _1 . This mass estimate is of 
the "enclosed" mass within the resolution of the interferometer, 
but since it is significantly higher than the best estimate of the 
envelope and disk masses of 0.026 and 0.0081 M respectively 
(see §|4]l, it is likely close to the total mass of the central star. 
This is also a lower limit since it is not corrected for inclination, 
but given the small minor axis relative to the major axis from 
the Gaussian fit to continuum emission, it is not unreasonable 
to assign a high inclination angle corresponding to an edge-on 
disk. 

In summary, we detect signatures of Keplerian rotation in the 
HCO + 3-2 emission for 4 Class I sources allowing direct deter- 
minations of their central stellar masses, which range from about 
0.3 to 2.5 M . As discussed in Appendix lA.2l the 100 AU scales 
of embedded protostars can only be traced with HCO + 3-2 in 



sources with envelope masses less than about 0.1 M , however. 
Future high sensitivity ALMA observations will be required to 
observe more optically thin isotopes, to determine the dynamical 
structure of the more deeply protostars. 

6. Discussion 

6.1. Effects of assumed dust properties and temperatures on 
disk masses 

As in most other studies in the literature, two assumptions in our 
study are that the disk can be characterized by a uniform tem- 
perature and by a single unchanging dust opacity law throughout 
its extent and evolution. Both these assumptions may affect the 
systematic trends observed in the data: for the temperature an 
increase in luminosity increases the mass-weighted temperature 
of the disk, while the shielding by the disk itself will lead to a 
lower mass weighted temperature with increasing mass. As an 
example, Fig. [14] compares the temporal evolution of the tem- 
perature at a radius of 200 AU between rays through the model 
with differen t inclinations with respect to the disk plane in the 
simulation of lVisser et all d2009h . Also shown is the evolution of 
the luminosity in the simulations. Early in the evolution, before 
the disk is formed, the temperature tracks the luminosity changes 
closely, but as the disk forms (at about 10 5 years) the tempera- 
ture at inclinations directed through the disk drops sharply, even 
though the luminosity continues to increase. The strong increase 
in total luminosity at about 3xl0 5 years is directly reflected in 
the temperature profiles. 

Fig. [T3] compares the relationships b etween (sub)m il limete r 
flux and disk ma ss adop ted by iTerebev et al.1 d 1993b . 
lAndrews & Williams! d2007bl) and |j0rgensen etail (|2007a) 
to the relationship predicted by the radiative trans fer modeling 
of the collapse simulation from lVisser et al.l (120091) : it is clearly 
seen that the models require a steeper relationship between the 
submillimeter flux and disk mass than what is obtained with a 
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Fig. 10. Archive 1.6 fim HST/NIC MOS images of TMR1, 
TMC1A, GSS30 and IRS 43 (see iTerebev et alJ 1 19981 and 
lAllen et al.l |2002| for details about the data) with HCO + 3- 
2 emission contours from the SMA overlaid. Panel a-c show 
TMR1, TMC1 A and GSS30 with blue and green contours corre- 
sponding to blue-shifted emission integrated from Vlsr -2AV to 
Vlsr - AV and V LSR - AV to V LSR with AV=1.0 km s _1 (TMR1 
and TMC1A) and AV=1.5 km s 1 (GSS30). For IRS 43 in panel 
d the blue and red contours correspond to blue- and red-shifted 
emission integrated from Vlsr-2AV to Vlsr - AV and Vlsr+AV 
to Vlsr + 2AV with AV=1.5 km s -1 . In that panel the white dot- 
ted lines indicate the rough extent of the dark lane in the HST 
image. 



single constant temperature: to calculate the disk mass from the 
submillimeter flux in the model, one would need a temperature 
scaling as: 



T = 11.3 K 



M di! 



0.0 IM, 



-0.25 



0.1 



(4) 



as shown with the solid line in Fig. [15]- i.e., an increasing tem- 
perature with increasing luminosity (and thus heating) and a de- 
creasing temperature with increasing disk mass (and thus shield- 
ing at a given radius). Although the variations in the temperature 
with disk mass and stellar luminosity are slow, they still intro- 
duce a systematic gradient in the derived disk masses by up to a 
factor of 4 through the evolution of the embedded stages of the 
YSOs. It thus appears that the disk masses in the later stages rel- 
ative to the early stages can be underestimated by about this fac- 
tor - solely due to the evolution of the protostellar system. This 
may explain the slightly higher masses for the Class sources in 
Fig. 13 

Fig.Q~6]compares the derived disk masses to bolometric tem- 
perature as well as the ratio of the source luminosity over the 
envelope mass, L*/M env . The latter ratio is an alternative evo- 
lutionary tracer, increasing with time as long as the YSO lumi- 
nosity is dominated by accretion and the stellar mass (and lumi- 
nosity) increases, whereas the envelope dissipates and its mass 
decreases. In this figure the Class disk masses are decreased by 





Fig. 12. Position-velocity (PV) diagram for the HCO + emission 
in IRS 43 directly from the image-plane (upper panel) and from 
fitting the individual channels in the (m, v)-plane (lower panel). 
The PV diagram has been measured along the major axis of the 
HCO + emission at a position angle of -70° (measured from north 
toward east). Overplotted are the predicted curves for Keplerian 
rotation around a 1.0 M Q central object (solid) line (best fit to 
(w, v)-data), which can be compared to the dashed lines showing 
the same predictions for 0.6 and 1.4 M Q central objects. 



a factor of 2 and the Class I disk masses increased by a factor of 

2 to simulate the temperature effect discussed above. With this 
correction, the Class and I disk masses are roughly constant 
with a median mass of about 0.04 M Q for the full sample. 

The assumption of a constant unchanging dust opacity law 
similar to the large-scale envelope dust may also be wrong: in 
particular, grain growth may increase the dust opacity at sub- 
millimeter wavelengths leading to a lower derived disk mass for 
a given interferometric flux. That something like this may be 
going on is hinted by the observation that the compact compo- 
nents around the Class sources show spectral slopes, a, of 
th e emission F v oc V, a t 230 and 345 GHz in the range 2- 

3 (|j0rgensen et al., 2007a). This suggests that the disks are ei- 
ther optically thick or have significantl y different dus t opacities 
than standard ISM dust. As discussed in J0r gensen et aU the disk 
masses implied for the Class objects suggest that only the most 
massive ones may be marginally optically thick. Resolved obser- 
vations at multiple wavelengths and measurements of their sizes 
are required to address this. On the other hand, the low spectral 
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Fig. 13. IRS 43: comparison between HCO + 3-2 spectrum from 
the SMA data integrated over the single-dish beam and JCMT 
observations. The interferometric spectrum has been scaled by a 
factor 2.5. 
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Fig. 14. Tempora l evolution of the t emperature and luminosity in 
the simulation of IVisser et al.l d2009l) . The solid lines shows the 
temperature at 200 AU along rays inclined with 0, 5, 10, 15, 20 
and 90° relative to the disk plane (from bottom to top; i.e., the 
line with the highest temperature corresponds to an inclination 
of 90° - perpendicular to the disk plane). The dashed line shows 
the corresponding evolution of the luminosity. 



indices confirm that the compact components around the Class 
sources are not associated with the larger scale envelopes for 
which an a of 3.5-4 is expected. 

Other dust opacity laws or t e mperat ures have been used in 
literature: lAndrews & Williamsl d2007bl) adopted an empirical 
scaling based on modeling of the disk SEDs, corresponding to 
an equiva lent dust temperature o f 20 K and a dust opacity-law 
following Bec kwith et a l. (1990) normalized to the dust opac- 
ity at 1000 GHz of 0.1 cm 2 g -1 (for a dust-to-gas ratio of 
1:100). The lower temperature leads to higher disk masses for 
a given 1.1 mm flux by a factor of 2.3 whereas their adopted 
dust opacity-law corresponds to lower disk masses b y a factor 
2.0 compared to ours. Likewise, Tereb ev et al.l (|T993) adopted a 
dust temperature of 40 K and adopted a dust opacity law with a 
mass opacity at 250 fim of 0.1 cm 2 g -1 . If compared at a com- 
mon wavelength of 850 /mi, these three choices of parameters 
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Fig. 15. Relationship between disk mass and 850 fim interfer- 
ometer flux assuming a distance of 250 pc: the three lower 

i 1 

straight lines indicat e the linear rel ationship s fromlTerebev et al.l 
d 19931) (lowest) lAndrews & Williamsl d2007bl) (middle) and 
Jol'ge nsen et alJ d2007ah (top) scaled to the same distance and 
wavelength. The filled symbols indicat e the relationship f ound 
for the radiative transfer modeling of the lVisser et al.l d2009l) col- 
lapse simulation discussed in the text. The solid line through 
these points indicates a relationship where the temperature used 
to relate the disk mass and flux follows Eq. [4] Finally, the three 
dotted lines indicate predicted line ar relationsh i ps usi ng the dust 
properties from the simulations of IVisser et al.l Q2009) with tem- 
peratures of 8, 11 and 16 K (from top to bottom). 



result in di sk masses that agree within + 15% from the middle 
value from lAndrews & Williamsl (l2007bl) . 

However, systematic changes in the dust properties occur- 
ring during the evolution of the young protostars could have 
an important effect on our derived disk masses and observed 
trends. For example, the dust opacities at 500 //m-l mm from 
lOssen kopf & Hennin g| (11994 increase by as much as an or- 
der of magnitude depending on the density at which the grains 
have coagulated, although the effect is significantly sma l ler for 
grains with ice mantles. As shown by ID'Alessio et alJ d200ll) 
grain growth increases the dust opacity at a specific wavelength, 
A, as long as the largest grains are smaller than A. In a sce- 
nario where gr ains grow toward size s inferred for T Tauri disks 
(~ 1mm; e.g.. ID'Alessio et al]|2001l) as the protostellar system 
evolve through its embedded stages the derived disk masses for 
the more evolved sources would be overestimated under our as- 
sumptions of an unchanging dust opacity law. This would work 
in the opposite direction of the temperature effect - toward even 
more massive disks in the Class compared to Class I stages. It 
cannot be ruled out that grain growth would appear in a fashion 
not straightforwardly increasing with "age" of the studied young 
stellar objects, e.g., if the observed dust in the Class sources is 
located in regions with relatively higher densities for a signif- 
icant fraction of time. Based on the current observational data, 
it is not possible to separate these scenarios, however: further 
resolved multi-wavelength images of the dust continuum emis- 
sion from the envelopes and disks around protostars or of well- 
calibrated gas tracers (taking their chemistry into account) are 
required. 
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Fig. 16. Disk masses as function of bolometric temperature f up- 
per) and source luminosity/envelope mass ratio (lower). In com- 
parison to Fig. [9] the Class disk masses have been decreased 
by a factor 2 and the Class I disk masses increased by a similar 
amount - simulating the evolutionary effect in the relationship 
between flux and disk mass suggested by Fig. [T31 



6.2. Collapse and disk formation models 

Power-law density profiles are simplified parameterizations of 
the underlying envelope density profile and the envelopes are 
not spherical due to, e.g., rotation and the action of protostellar 
outflows. 

In the framework of the Isiiuld 19771) and lTerebev et al.1 dl984 
models for inside-out collapsing envelopes, the density profile 
flattens from a stationary r~ 2 density profile through a free- 
falling r~' 5 profile to a flattened inner en velope where t he an- 
gular momentum roughly balances gravity. Terebev etld] d 1993b 
examined the millimeter emission from such envelope models 
using analytical descriptions for the envelope density and tem- 
perature profiles and used those to compare IRAM 30 m single- 
dish continuum observations (11"; at 1.3 mm) to OVRO in- 
terferometric observations (4-10"; at 2.7 mm). Consistent with 
Fig.|6]they found that the emission in the interferometric obser- 
vations is completely dominated by a central compact compo- 
nent whereas most of the observed single-dish flux in fact has 
its origin in the collapsing envelope. Th e dominance o f the disk 
to the interferometric flux was found by iTerebev et alj to be due 
to ( i) the smaller beam of the interferometric data and negligible 
envelope mass on those scales and (ii) the steeper wavelength 



dependence found for the dust emission in the optically thin en- 
velope. 

We expand this analysis by introducing the formation of 
a circ umstellar vi s cous disk in the center using the mod- 
els of IVisser et al.l d2009l) and following the evolution of the 
subm illimeter b rightness profiles observed by an interferome- 
ter: Visser ^tail adopted t he framework of ICasse n & Moosman 
( 1981) and ITerebev et al.l d 19841) to describe a collapsing enve- 
lope and introduced the formation and gr owth of a 2-dimensional 
centra l disk and star similar to, e.g., lYorke & Bodenheimerl 
( 1999). The models of Visse r et ail also include the presence of 
an outflow cavity with a growing openin g angle through the evo - 
lution consistent with observations (e.g jArce & Sargenl l2006). 
Using the axisymmetric radi a tive transfer code "RADMC" 
dDullemond & Dominikl 120041) . I Visser et al.l calculated the tem- 
perature profile for time-steps in the above model description 
and furthermore followed the chemical history of particles dur- 
ing the collapse and disk formation. For this purpose, we adopt 
the temperature and densi ty profiles o f the envelope (with out- 
flow cavity) and disk from lVisser et al.1 as well as their prescrip- 
tion for the evolution of the stellar luminosity and effective tem- 
perature and dust parameters, and ray-trace these. Fig.fTTlshows 
the evolution of the predicted visibility curves at di fferent time- 
steps for two representative models from the grid of lVisser et al.l 

Both panels clearly illustrate the effect discussed above: un- 
til the formation of the circumstellar disk more than 95% of the 
total flux at a given wavelength is emitted on baselines shorter 
than 50 kA - corresponding to a ratio between the 1 . 1 mm inter- 
ferometric flux and the peak 850 pm single-dish flux of less than 
a few % (depending on the exact spectral slope between the ob- 
servations at these differing wavelengths) in accordance with the 
results above. Once the disk is formed, the longer baseline flux 
rapidly increases - first simultaneously on all baselines and then 
more rapidly on the shorter baselines as the disk grows to sizes 
large enough to be resolved by the interferometric observations. 
It should be noted that even though the outflow cavity is opened 
up early in the evolution of the collapsing protostar, it does not 
provide significant compact emission, e.g., due to heated ma- 
terial in the outflow cavity walls, to affect the evolution of the 
observed interferometric flux. 

An alternative explanation to the disk hypothesis favored 
in this work is that there is an enhancement of the den- 
sity in the envelop e profile at small scales. As an example, 
Chiang et al. (2008) ad opted the theoretical collapse models of 
Tassis & Mouschovias (2005 ) and compared t hose to 2.7 mm 
interferometric data from [Lqonev et alj d2000l) . The models of 
Tassis & Mous chovias! d2005l) simulate the formation of a pro- 
tostellar object from a cloud core in the presence of magnetic 
fields. In their simulation, ambipolar diffusion leads to the for- 
mation of a thermally and magnetically super-critical core which 
contracts dynamically. Initially this core has a density profile 
n oc r~ lJ , but as mass and magnetic flux is accreted from the 
envelope onto the central protostar, a region of enhanced mag- 
netic fields is formed driving an outward-propagating shock and 
causing a density enhancement growing quasi-periodically from 
scales of ~100 AU. This density enhancement on small scales 
will result in a compact component observable by the interfer- 
ometer which is slightly resolved on long baselines and thereby 
in some cases allevi ates the need for a ci rcumstellar disk on those 
scales as shown by Chian get al.l ((2008). 

Still, the conclusions of both studies are similar: a density en- 
hancement is required in the collapsing envelope on small scales 
- above what would be expected from, e.g., a free-falling enve- 
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Fig. 17. Evolution of the predicted visibility amplitude curves 
(normalized by the zero-spacing flux) as function of projected 
baseline lengths for two models fo r envelope collapse and disk 
formation from IVisser et"aT] d2009h characterized by the solid- 
body rotation rate, Qo, and the effective sound speed c s . Upper: 
O = 1(T 14 s _1 and c s = 0.19 km s -1 . Lower: Q = 10~ 13 s _1 
and c s - 0.26 km s . The different curves correspond to snap- 
shots at different times indicated in the upper right insert in each 
panel, also showing the evolution of the envelope, disk and stel- 
lar masses. 



lope. The question is whether this component is inherent to the 
collapse profile or the first sign of the formation of a circumstel- 
lar disk. A straightforward, although currently observationally 
challenging, way to address this issue is through resolved line 
observations of the innermost envelope: the resolved line obser- 
vations of the Class I sources in this sample dBrinch et alJl 2007; 
lLommen et all 120081) show clear signatures of Keplerian rota- 
tion, which argues in favor of the disk scenario for the sources in 
those evolutionary stages at least. For the more deeply embed- 
ded protostars future ALMA observations are needed to confirm 
the dynamical nature of the central compact dust component. 



6.3. Constraints on evolutionary models 

The presented observations and methodology make it possible 
to directly constrain models for the evolution of young stellar 
objects - in particular, concerning the transport of matter from 
the larger scale envelope to the central star. This is in principle 
also the most direct way for estimating the evolutionary status 
of individual young stellar objects - rather than indirect indica- 
tors such as their spectral energy distributions or similar. Fig. ITSI 
compares the predicted stellar and disk masses relative to the 
envel ope masses ( i .e., M star /M env vs. M^/M em ) from the mod- 
els of IVisser et al.l (120091) to our observed values. The thickened 
part of the lines in this plot indicates where the Class and I 
sources would fall on the evolutionary track given their mea- 
sured Mdi s k/M env ratios. 

For the observed ratios of M^IM &W for the Class sources, 
the models predict ratios of stellar to envelope masses M star /M env 
in the range « 0.2 - 0.6, i.e., systems that have accreted « 15- 
40% of their final mass - in accordance with the phenomenologi- 
cal definition of objects in the Class stage. Likewise the Class I 
sources are predicted to have more massive stars than envelopes 
- again, in accordance with the phenomenological distinction 
between Class and I objects. For the Class I sources with mea- 
sured stellar masses, the observed disk-to-stellar mass ratios are 
in the 1-10% range, whereas the models predict a higher ratio of 
20-25% in the late Class I stage (i.e., overestimate the disk mass 
relative to the stellar mass). It is seen that the predicted stellar 
masses for the Class sources with this model are a fewxO. 1 M Q . 

The low disk masses relative to stellar masses suggest that 
material is quickly and efficiently accreted from the envelope 
through the disk onto the star. This is also directly illustrated 
by the comparison between the observed and modeled stellar- 
to-disk mass ratios for the five Class I sources for which both 
have been determined. For the Class I sources the observed disk 
masses are 1-10% of the corresponding stellar masses whereas 
the models predict ratios of 20-25%. This suggests that the ma- 
terial is processed more quickly t hrough the circums tellar disks 
than what is given by the models. Loon ev et alJ d2003) suggested 
that this was the case based on their conclusion that Class pro- 
tostars have disks that are not more massive than their Class I 
counterparts. Our results are in agreement with this general con- 
clusion. 

It is worth emphasizing though, that our analysis does 
find evidence for disks around the Class sources. If the 
Class sources indeed are precursors to the Class I sources 
in this sample, this in turn would suggest that circumstellar 
disks are formed early in the evolution of low-mass proto- 
stars. Without resolved line observations it is not possible to 
address whether these disks in fact are rotationally supported 
Keplerian disks. The absence of a c lear rotation signat ure in 
the disk around NGC 1333-IRAS2A dBrinch et alll2009l) com- 
bined with its large apparent size fr om the dust continuum mea- 
surements (J0rgens en et all 12005 af) could perhaps be consistent 
with it being an uns table "pseudo-disk" as in the framework of 
lGalli&Shuldl993allbh - although it may also simply be that the 
evolutionary time-scale for IRAS2A sti ll has been too sho rt for 
the disk to establish Keplerian rotation (Bri nch et all 12009). 

An open question in the comparison to the theoretical mod- 
els is how important outflows are in carrying away envelope 
material. The observed increasing opening angles of outflows 
with increasing age of the protostar has been interpreted as 
evidence for the erosion of the protostellar envelope b y the 
outflow dVelusamv & LangeA 119981: lArce & Sargentl [2006h in 
part explaining the observed shift in the peak between the core 
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Fig. 18. Predicted stellar mass, M star , vs. disk mass, M^, both 
measured relative to the envelope mass M env in the models of 
IVisser et al.ld2009l) with Q» = 10~ 14 s _I and c s of (a) 0.19 km s" 1 
and (b) 0.26 km s _! (solid lines; see also Fig.fTTll. On this model 
track, the thickened lines indicate where the Class and I objects 
would be located given their distribution of measured M^/M^ 
values (red and blue, respectively). The sources for which stellar, 
disk and envelope masses are measured are shown with symbols. 
Finally, the dashed lines indicate Mdisk/A^star ratios of 1 % (upper) 
and 10% (lower). 



and initial mass fu nction (e.g., iMversl 120081) . In the models of 
Viss er et af] ((2009) and shown in Fig. [18] the amount of ma- 
terial carried away by the outflow only amounts to a few per- 
cent. Generally, a more efficient outflow dispersal of the core 
would decrease the disk and stellar masses relative to the enve- 
lope masses - corresponding to a diagonal shift of a given mod- 
eled object in Fig. [18] - roughly along the evolutionary tracks. 
Therefore, unless the outflow preferentially alter the accretion 
onto the disk relative to the star, the ratio between the masses of 
the two will not be changed and thus also not the implications 
discussed here. 

A final interesting possibility is that the flat distribution of 
disk masses from the Class to I stages reflects the stability 
of the disks, given the suggestion that the luminosity problem 
of low-mass protostars can be explained if mass is accreted on 
large scales in the disks and from there onto the central star in 
highly episodic events: in their models for c ollapse of a dense 
core a nd formation of a protostar and disk, Vorobvov & Basu 
(2006) find strong episodic accretion once the disk has formed 
which keeps its mass lower than the stellar mass at all points dur- 
ing the protostellar evolution. Given that the recent Spitzer/c2d 
results suggest that the luminosity pro blem extends to the earli- 
est Class stages dEvans et all l2009h . it is tempting to suggest 
that material already at these stages indeed ends up in circum- 
stellar disks and from there is accreted toward the central star in 
bursts when a local instability criterion is exceeded - e.g., related 
to comparison between the disk size/mass and stellar mass. This 
could in turn flatten the distribution of disk masses as function of 
bolometric temperatures, with an upper envelope of disk masses 
of 0.05-0. 1M corresponding to the maximum values observed 
for our sample of Class and I protostars. It should be empha- 
sized that since our sample is rather small, biases may have been 
introduced in its selection. Future studies of larger samples of 
protostars with current and upcoming (sub)millimeter interfer- 
ometers will be needed to address this potential bias. 



7. Summary and outlook 

We have presented a comprehensive analysis of the dust contin- 
uum envelope emission of a sample of 20 embedded low-mass 
YSOs (Class and I) based on submillimeter interferometric 
observations from the Submillimeter Array probing scales of a 
few hundred AU and single-dish continuum observations from 
the JCMT (scales larger than 2000 AU). In addition we discuss 
HCO + 3-2 line emission for the Class I sources in the sample. 
With these data we have characterized the evolution of the three 
main components YSOs - the envelope, disk and stellar masses 
- and the relation between these compared to theoretical models. 
The main conclusions are as follows: 

1 . The millimeter dust continuum emission from both the in- 
terferometric observations (on scales of 250-500 AU) and 
single-dish observations (scales larger than 2000 AU) de- 
creases in strength from the more deeply embedded Class 
to the more evolved Class I sources. Relative to the single- 
dish data, the interferometric continuum flux at millimeter 
wavelengths increases toward the more evolved sources. 

2. Through comparison to predictions from dust radiative trans- 
fer models, an upper limit on the contribution from proto- 
stellar envelope to the interferometric continuum flux is es- 
timated - relative to the continuum flux measured on larger 
scales with single-dish telescopes. For the Class I sources 
about 87% of the observed SMA continuum flux (median 
value for the sample of Class I sources) has its origin in 
the circumstellar disks with the remainder coming from the 
large-scale envelope. For the Class sources, 68% of the in- 
terferometric emission is related to the circumstellar disks. 

3. Using empirical relationships between the envelope and disk 
fluxes and their masses from the radiative transfer models, 
the observations show that the disk-to-envelope mass ratio 
increases from about 1-10% in the Class stage to about 
20-60% in the Class I stage, reflecting the dissipation of the 
envelope in the more evolved stages. 

4. The compact interferometric fluxes due to disks are observed 
to decrease from the Class to the Class I objects. However, 
this decrease is largely due to the changing physical prop- 
erties, in particular disk temperature, during this phase. If 
these effects are taken into account, the disk mass shows no 
significant trends from the Class to I stages and is typically 
0.05 Mq, albeit with significant scatter. Thus, no evidence is 
seen for significant disk mass growth from the Class to the 
Class I stages. 

5. Four Class I sources in the sample show evidence of 
Keplerian rotation in HCO + 3-2, and for IRS 43 also in HCN 
3-2. IRS 63 and IRS 43 show evidence for elongated "con- 
tinuum disks" at 1.1 mm. For IRS 43 both the "continuum 
disk" and the HCO + 3-2 emission coincides with a dark lane 
in the HST 1.6 yum images - and orientated perpendicular to 
the direction of the protostellar outflow. For the four Class 
I sources we derive central masses of 0.3-2.5 M Q from the 
Keplerian rotation curve. 

6. The high excitation HCO + and HCN lines appear to be good 
tracers of the motions in the inner envelopes and disks, al- 
though in the most massive envelopes their use for this pur- 
pose is limited by the lines of the main isotopes being opti- 
cal thick. Future observations with ALMA will allow similar 
studies using the less abundant isotopologues. 

7. We compare the derived envelope and stellar masses to mod- 
els for the mass evolution in a collapsing protostellar sys- 
tem includi ng the format i on and growth of a central circum- 
stellar disk dVisser et al.L 120091) . The data are in agreement 
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with phenomenological definition of Class sources having 
accreted less, and Class I sources more, than half of their 
final mass. The models cannot explain the observed stellar 
masses for the handful of Class I sources for which these 
have been measured from the disk dynamical structure: gen- 
erally the models produce too massive disks relative to their 
stellar masses. 

8. Taken together these results suggest that circumstellar disks 
are established early in the evolution of low-mass proto- 
stars, but at the same time transport material rapidly from the 
larger scale envelopes onto the central stars. This may occur 
in episodic events, which may explain why the "luminosity 
problem" (that embedded YSOs are under-luminous com- 
pared to the theoretical predictions from pure accretion mod- 
els) a ppears already du ring the Class protostellar stages 
(e.g JEvans et al.Ll2009h . 

An important test of the theoretical models for the evolution 
of embedded young stellar objects will be to measure the disk 
dynamical structure even in the most deeply embedded stages. In 
this way it will be possible to address whether these disks have 
in fact settled into a rotationally stable pattern - and if so what 
central stellar masses are for the central newly formed star - or 
if the disk dynamics perhaps suggest some other configuration, 
for example, unstable "pseudo-disks". In this context, it remains 
a fundamental question whether the central stellar masses of the 
Class objects indeed are significantly lower than those of the 
Class I sources, which would confirm their earlier evolutionary 
status. 

Future ALMA observations will provide the sensitivity to 
molecules that are less abundant, and thus not severely opti- 
cally thick in the outer envelopes, and thus a tracer of this struc- 
ture. With ALMA it will also be possible to perform larger sur- 
veys of a statistically more unbiased sample of YSOs to con- 
strain the evolution more directly. This will for example address 
whether traditional classification schemes for young stellar ob- 
jects need to be recalibrated, whether the episodic accretion sce- 
narios for YSOs are tied to the properties of their circumstellar 
disk, whether there are significant variations between accretion 
histories of protostars in different star forming environments, 
and what fraction of the envelope material is accreted onto the 
star and what fraction is carried away by the outflows. 

Acknowledgements. We thank the referee for good suggestions about the pre- 
sentation of these results. It is also a pleasure to thank the Submillimeter Array 
staff for help in carrying out observations and technical assistance with the re- 
sulting data. The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of Astronomy and 
Astrophysics, and is funded by the Smithsonian Institution and the Academia 
Sinica. This research was based in pail on observations from the James Maxwell 
Telescope archive at the Canadian Astronomy Data Centre operated by the 
National Research Council of Canada with the support of the Canadian Space 
Agency. It is also used observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope Institute. STScI 
is operated by the association of Universities for Research in Astronomy, Inc. 
under the NASA contract NAS 5-26555. Astrochemistry in Leiden is supported 
by a Spinoza Grant from the Netherlands Organization for Scientific Research 
(NWO) and a NOVA grant. The research of TLB was supported by NASA 
Origins grant NNX09AB89G. 



References 

Allen, L. E., Myers, P. C, Di Francesco, J., et al. 2002, ApJ, 566, 993 

Andre, P. & Montmerle, T. 1994, ApJ, 420, 837 

Andre, P., Ward-Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 

Andre, P., Ward-Thompson, D., & Barsony, M. 2000, in Protostars and Planets 

IV, ed. V. Mannings, A. P. Boss, & S. S. Russell (University of Arizona Press, 

Tucson), 59 



Andrews, S. M. & Williams, J. P. 2005, ApJ, 631, 1 134 
Andrews, S. M. & Williams, J. P. 2007a, ApJ, 671, 1800 
Andrews, S. M. & Williams, J. P. 2007b, ApJ, 659, 705 
Arce, H. G. & Sargent, A. I. 2006, ApJ, 646, 1070 

Barsony, M., Kenyon, S. J., Lada, E. A., & Teuben, P. J. 1997, ApJS, 112, 109 
Basu, S. 1998, ApJ, 509, 229 

Beckwith, S. V. W., Sargent, A. I., Chini, R. S., & Guesten, R. 1990, AJ, 99, 924 
Brinch, C, Crapsi, A., J0rgensen, J. K., Hogerheijde, M. R., & Hill, T. 2007, 
A&A, 475, 915 

Brinch, C, J0rgensen, J. K., & Hogerheijde, M. R. 2009, A&A, 502, 199 
Brown, D. W. & Chandler, C. J. 1999, MNRAS, 303, 855 
Cassen, P. & Moosman, A. 1981, Icarus, 48, 353 

Chandler, C. J., Terebey, S., Barsony, M., Moore, T. J. T., & Gautier, T. N. 1996, 
ApJ, 471,308 

Chiang, H.-E, Looney, L. W., Tassis, K., Mundy, L. G., & Mouschovias, T. C. 

2008, ApJ, 680, 474 
Curiel, S., Girart, J. M., Rodriguez, L. E, & Canto, J. 2003, ApJ, 582, L109 
D'Alessio, P., Calvet, N., & Hartmann, L. 2001, ApJ, 553, 321 
Di Francesco, J., Johnstone, D., Kirk, H, MacKenzie, T, & Ledwosinska, E. 

2008, ApJS, 175, 277 
Duchene, G., Bouvier, J., Bontemps, S., Andre, P., & Motte, F. 2004, A&A, 427, 

651 

Dullemond, C. P. & Dominik, C. 2004, A&A, 417, 159 
Dunham, M. M., Crapsi, A., Evans, II, N. J., et al. 2008, ApJS, 179, 249 
Eisner, J. A., Hillenbrand, L. A., Carpenter, J. M., & Wolf, S. 2005, ApJ, 635, 
396 

Enoch, M. L., Glenn, J., Evans, II, N. J., et al. 2007, ApJ, 666, 982 
Evans, N. J. 1999, ARA&A, 37, 311 

Evans, N. J., Dunham, M. M., J0rgensen, J. K., et al. 2009, ApJS, 181, 321 
Galli, D. & Shu, F. H. 1993a, ApJ, 417, 220 
Galli, D. & Shu, F. H. 1993b, ApJ, 417, 243 

Girart, J. M., Curiel, S., Rodriguez, L. E, et al. 2004, AJ, 127, 2969 
Girart, J. M., Rao, R., & Estalella, R. 2009, ApJ, 694, 56 
Girart, J. M., Rodriguez, L. E, & Curiel, S. 2000, ApJ, 544, L153 
Grosso, N, Alves, J., Neuhauser, R., & Montmerle, T. 2001, A&A, 380, LI 
Gutermuth, R. A., Myers, P. C, Megeath, S. T, et al. 2008, ApJ, 674, 336 
Haisch, Jr., K. E., Greene, T. P., Barsony, M., & Stahler, S. W. 2004, AJ, 127, 
1747 

Hatchell, J. & Fuller, G. A. 2008, A&A, 482, 855 

Hatchell, J., Fuller, G. A., Richer, J. S., Harries, T. J., & Ladd, E. F. 2007, A&A, 
468, 1009 

Ho, P. T. P., Moran, J. M., & Lo, K. Y. 2004, ApJ, 616, LI 
Hogerheijde, M. R. & Sandell, G. . 2000, ApJ, 534, 880 

Hogerheijde, M. R., van Dishoeck, E. E, Blake, G. A., & van Langevelde, H. J. 

1997a, ApJ, 489, 293 
Hogerheijde, M. R., van Dishoeck, E. E, Blake, G. A., & van Langevelde, H. J. 

1998, ApJ, 502, 315 

Hogerheijde, M. R., van Langevelde, H. J., Mundy, L. G., Blake, G. A., & van 

Dishoeck, E. F. 1997b, ApJ, 490, L99+ 
Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, User Manual for DUSTY, 

University of Kentucky Internal Report 
J0rgensen, J. K. 2004, A&A, 424, 589 

J0rgensen, J. K, Bourke, T. L., Myers, P. C, et al. 2007a, ApJ, 659, 479 
J0rgensen, J. K, Bourke, T. L., Myers, P. C, et al. 2005a, ApJ, 632, 973 
J0rgensen, J. K, Harvey, P. M., Evans, N. J. I., et al. 2006, ApJ, 645, 1246 
J0rgensen, J. K., Hogerheijde, M. R., van Dishoeck, E. E, Blake, G. A., & 

Schoier, F. L. 2004a, A&A, 413, 993 
J0rgensen, J. K, Johnstone, D., Kirk, H, & Myers, P. C. 2007b, ApJ, 656, 293 
J0rgensen, J. K, Johnstone, D., Kirk, H, et al. 2008, ApJ, 683, 822 
J0rgensen, J. K, Lahuis, F, Schoier, F. L., et al. 2005b, ApJ, 631, L77 
J0rgensen, J. K, Schoier, F. L., & van Dishoeck, E. F. 2002, A&A, 389, 908 
J0rgensen, J. K, Schoier, F. L., & van Dishoeck, E. F. 2004b, A&A, 416, 603 
Keene, J. & Masson, C. R. 1990, ApJ, 355, 635 

Kenyon, S. J., Hartmann, L. W., Strom, K. M., & Strom, S. E. 1990, AJ, 99, 869 
Lay, O. P., Carlstrom, J. E., Hills, R. E., & Phillips, T. G. 1994, ApJ, 434, L75 
Leous, J. A., Feigelson, E. D., Andre, P., & Montmerle, T. 1991, ApJ, 379, 683 
Lommen, D., J0rgensen, J. K., van Dishoeck, E. F, & Crapsi, A. 2008, A&A, 
481, 141 

Looney, L. W., Mundy, L. G., & Welch, W. J. 2000, ApJ, 529, 477 
Looney, L. W., Mundy, L. G., & Welch, W. J. 2003, ApJ, 592, 255 
Looney, L. W., Tobin, J. J., & Kwon, W. 2007, ApJ, 670, L131 
Montmerle, T., Grosso, N, Tsuboi, Y, & Koyama, K. 2000, ApJ, 532, 1097 
Mundy, L. G, Looney, L. W., Erickson, W., et al. 1996, ApJ, 464, L169 
Myers, P. C. 2008, ApJ, 687, 340 
Myers, P. C. & Ladd, E. F. 1993, ApJ, 413, L47 
Ossenkopf, V. & Henning, T. 1994, A&A, 291, 943 

Padgett, D. L., Rebull, L. M., Stapelfeldt, K. R., et al. 2008, ApJ, 672, 1013 
Qi, C. 2008, The MIR Cookbook, The Submillimeter 



20 



J. K. J0rgensen et al.: The mass evolution of embedded protostars 



Array / Harvard-Smithsonian Center for Astrophysics 

( http://cfa-www.hai'vard.edu/~cqi/mircook.htmli 
Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Data 

Analysis Software and Systems IV, ed. RA. Shaw, H.E. Payne and JJ.E. 

Hayes, PASP Conf Series 77, 433 
Schoier, F. L., J0rgensen, J. K., van Dishoeck, E. F., & Blake, G. A. 2002, A&A, 

390, 1001 

Schoier, F. L., J0rgensen, J. K., van Dishoeck, E. F, & Blake, G. A. 2004, A&A, 
418, 185 

Shirley, Y. L., Evans, N. J., & Rawlings, J. M. C. 2002, ApJ, 575, 337 
Shu, F. H. 1977, ApJ, 214, 488 

Stutz, A. M., Rubin, M., Werner, M. W., et al. 2008, ApJ, 687, 389 
Tassis, K. & Mouschovias, T. C. 2005, ApJ, 618, 783 
Terebey, S., Chandler, C. J., & Andre, P. 1993, ApJ, 414, 759 
Terebey, S., Shu, F. H, & Cassen, P. 1984, ApJ, 286, 529 
Terebey, S., Van Buren, D., Matthews, K., & Padgett, D. L. 2000, AJ, 1 19, 2341 
Terebey, S., van Buren, D., Padgett, D. L., Hancock, T, & Brundage, M. 1998, 
ApJ, 507, L71 

Tobin, J. J., Hartmann, L., Calvet, N., & D'Alessio, P. 2008, ApJ, 679, 1364 
Tsuboi, Y., Imanishi, K., Koyama, K., Grosso, N., & Montmerle, T. 2000, ApJ, 

532, 1089 
Ulrich, R. K. 1976, ApJ, 210, 377 

van der Tak, F. F. S., Black, J. H., Schoier, F. L., Jansen, D. J., & van Dishoeck, 

E. F. 2007, A&A, 468, 627 
van Kempen, T. A., van Dishoeck, E. F, Salter, D. M., et al. 2009, A&A, 498, 

167 

Velusamy, T. & Langer, W. D. 1998, Nature, 392, 685 

Visser, R., van Dishoeck, E. F, Doty. S. D., & Dullemond, C. P. 2009, A&A, 
495, 881 

Vorobyov, E. I. & Basu, S. 2006, ApJ, 650, 956 

Ward-Thompson, D. 2007, Joint 32nd International Conference on Infrared 
and Millimeter Waves, and 15th International Conference on Terahertz 
Electronics. IRMMW-THz (ed. M. Griffin, P. C. Hargrave, T. J. Parker and 
K. P. Wood)., 720 

Weintraub, D. A., Kastner, J. H, Griffith, L. L., & Campins, H. 1993, AJ, 105, 
271 

Yorke, H. W. & Bodenheimer, P. 1999, ApJ, 525, 330 

Young, C. H, Shirley, Y. L., Evans, N. J., & Rawlings, J. M. C. 2003, ApJS, 145, 
111 

Young, E. T, Lada, C. J., & Wilking, B. A. 1986, ApJ, 304, L45 
Young, K. E., Enoch, M. L., Evans, II, N. J., et al. 2006, ApJ, 644, 326 
Zhang, Q., Wootten, A., & Ho, P. T. P. 1997, ApJ, 475, 713 



J. K. J0rgensen et al.: The mass evolution of embedded protostars, Online Material p 1 



Online Material 



J. K. J0rgensen et al.: The mass evolution of embedded protostars, Online Material p 2 



Appendix A: The use of HCO + as a dense gas tracer 

A.1. Outflow cavities in TMR1, TMC1A and GSS30-IRS1 

The HCO + 3-2 alignment with the near-infrared scattering nebu- 
losities in the HST images for TMR1, TMC1A and GSS30-IRS1 
(Fig.fTUb is similar to what has been found in high angular reso- 
lution o bservations of the lower ex citation 1-0 transition in both 
Class I dHogerheiide et al.Lll997at) and Class dJ0rgensenll 2OO4) 
sources. 

The patterns in the velocity maps (Fig.fTTI) support the inter- 
pretation that the HCO + 3-2 emission has its origin in one-side 
of the outflow cone: in all maps the HCO + emission is blue- 
shifted with respect to the systemic velocity with the most ex- 
treme velocities (most blue-shifted) narrowly confined around 
the continuum position and material closer to the systemic veloc- 
ity extending around this. This signature can best be understood 
in a simple geometry where HCO + probes material in a large in- 
falling protostellar envelope being swept up by the outflow, as 
illustrated in Fig. IA.1I If the velocities in the swept-up material 
are relatively small (comparable to the infalling velocities in the 
envelope) and HCO + 3-2 optically thick, the front-side of the 
envelope (toward us) will be slightly red-shifted due to the in- 
fall and thus obscure the outflowing material in the red-shifted 
cone on the far-side of the envelope (Fig. IA.U . which in turn 
is resolved-out. In contrast, the outflow cone pointed toward us 
will be unobscured by the envelope. Likewise, the near-infrared 
scattered light from the red-lobe would largely be blocked - in 
contrast to the scattered emission from the outflow cavity on the 
front-side. An alternative explanation is that all the sources are 
at the far-side of their respective clouds, the red-shifted part of 
the outflow cone breaking out into a low-density medium not 
observable in the line tracers. This would require a rather sharp 
density gradient, however. 

A.2. Limitations ofHCO + 3-2 as a tracer of disks 

The results above suggest that HCO + 3-2 is a good tracer of mo- 
tions in disks for envelopes that have already been largely dis- 
persed. For the more embedded YSOs, HCO + 3-2 emission may 
be of limited use as a tracer of rotation in the disk due to its opti- 
cal thickness. A simple estimate can be made of the typical enve- 
lope mass below which HCO + 3-2 line emission is optically thin 
in the envelope and can be used as a tracer of the kinematics in 
the innermost regions of the protostar. First, for typical envelope 
parameters we expect temperatures in the range of 15^-0 K and 
densi ties of 10 4 -10 7 cm" 3 . Us ing the Radex escape probability 
code ( Ivan der Tak et all |2007|PI for these parameters, the HCO + 
3-2 line has an optical thickness « 1 for an HCO + column den- 
sity of 5 x 10 12 cm" 2 and a line width of 1 km s . Secondly, we 
can estimate what density, hiooau is required to reach that opti- 
cal thickness in a free-falling envelope with a power-law density 
profile, n - (r/100AU)" 15 at a size scale of 100 AU correspond- 
ing to the region where the disk potentially forms, i.e., by simply 
estimating «iooau from 

/~10 4 AU 

N H co + = mooAu(r/100AU)- L5 dr (A.l) 

Jl0 2 AU 

with A^hco+ = 5 x 10 12 cm" 2 . Assuming a ty pical abundance 
of HCO + of 1 x 10" 9 with respect to H2 (e.g.. J0r gensen et"ail 
2004b), we find hiooau = 2 x 10 6 cm" 3 , which in turn corre- 
sponds to an envelope mass of 0.1 M Q . This is the maximum 

3 http : //www . sron . rug . nl/ ~vdtak/ radex/radex . php 
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Fig. A.l. Schematic figure of a collapsing envelope with an out- 
flow sweeping up material. The dashed lines indicate the r — 1 
surfaces for a red- and blue-shifted velocity (r rec i and Tt,i ue , re- 
spectively) for a line which becomes optically thick at low col- 
umn densities through the envelope. As illustrated in the figure, 
the blue-shifted part of the outflow cone may appear in front of 
the Tbiue = 1 surface, while the red-shifted part is behind the 
Tied = 1 surface, obscured by material in the large scale collaps- 
ing envelope, which in turn may be resolved out by the interfer- 
ometer. 

envelope mass for which we can use the HCO + 3-2 to trace 
scales smaller than about 100 AU without worrying about opac- 
ity effects in the larger scale envelope. It indeed corresponds to 
the estimate of the masses of the most deeply embedded Class I 
sources in our sample for which the outflow emission is seen. 
For these sources and the even more deeply embedded Class 
sources a less abundant or more optically thin tracer is therefore 
required to trace the dynamical structure on less than 100 AU 
scales. Future observations with ALMA will provide the sensi- 
tivity to image, e.g., the optically thin 3-2 and 4-3 submillimeter 
transitions of H 13 CO + in the disks around even more deeply em- 
bedded protostars and thus provide stellar masses for protostars 
in the deeply embedded stages. 

Appendix B: Summary of observed Class I sources 

B.1. IRS 43 

IRS 43 (YLW15; lYoung et all [l986h has been the target of 
numerous studies across the full wavelength range. At radio 
wavelengths, IRS 43 is resolved into two separate compo- 
nents, YLW15-VLA1 and -VLA2 with a separation of 0.6" 
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dGirart et all 120001) . From multi-epoch high angular resolution 
radio data, these binary components were found to show relative 
proper motion s indicative of or b ital motions in a 1.7 M total 
mass binary (Curie l et al.L 120031) . iGirart et alj found that one of 
the two components, VLA1, were resolved in observations at 3.6 
and 6 cm with the VLA with a deconvolved size of about 0.4" 
along its major axis and suggested that it was associated with a 
thermal radio jet. In addition, IRS 4 3 is a peculiar X-ray emitter 
showing quasi-periodic X-ray flares (Tsub oiet al.Ll2000r) consis- 
tent with magnetic shearing and reconnection between the cen- 
tral star and an accretion disk, with the period su ggesting a mass 
of the central star in the range 1.8-2.2 M Q dMontmerle et al.L 
l2000h . 

Th e region was imaged in the near-infrared bv lGrosso et alj 
(2001) who discovered a set of "embedded" Herbig-Haro 
objects, well aligned with the thermal jet candid ate from 
the centimeter observations of Gir art et ail d2000l) . Spitzer 
Space Telescope images of the embedded Herbig-Haro objects 
(Fig. IB. II ) show the knots clearly along a line pointing directly 
back toward the group of objects where IRS 43 and another 
embedded pr otostar, IRS 44, is located. As also pointed out by 
iGrosso et all the association of the Herbig-Haro knots to IRS 43 
is not unique: the nearby protostar, IRS 54 (see also jjB.5l ). 
also drives a large-scale, precessing outflow clearly seen in the 
Spitzer images in Fig. IB. II 

The dotted lines in Fig. IB. 11 show the propagation of the 
IRS 54 outflow, following the H2 emission in the eastern lobe 
and assuming that the precession of the IRS 54 outflow is sym- 
metric around the protostar - thereby giving the expected loca- 
tion of the western lobe. The western lobe in this case passes 
through the location of the southern part of the near-infrared 
Herbig-Haro objects: based on the ir morphologies in the higher 
resolution near-infrared images, IGrosso et alj d200ll) suggested 
that these knots (crosses in the lower left panel of Fig. IB.lt 
are indeed caused by the IRS 54 outflow with the northern 
knots (plus signs in the lower left panel of Fig. IB. II related to 
the IRS 43 outflow - the lower spatial resolution Spitzer im- 
ages do not reveal this distin ction clearly, though. Some sup- 
port for the interpretation of IGrosso et alj can be taken from 
the morphology of another bow-shock seen in the Spitzer im- 
ages at (a,(5)j2000 = (16 h 27 m 39.5 s , -24°33'02") to the west of 
the near-infrared Herbig-Haro knots. Together with the near- 
infrared Herbig Haro knots detected by IGrosso et all this bow- 
shock is on the southern side of the dashed-line suggesting that 
the dashed line indeed is delineating the northern edge of the 
western lobe (and the southern edge of the eastern lobe) of the 
IRS 54 outflow. It is also possible that another source in the 
group encompassing IRS 43 could be responsible for the Herbig- 
Haro knots: the main argument for the association to IRS 43 is 
the coincidence with the radio jet and the identification of a bow- 
shock south of IRS 43 in the near-infrared, possibly represe nting 
the counter-jet dGrosso et all 120011: iGirartetall l20odl2004l) . 

We detect IRS 43 in both HCO + 3-2, HCN 3-2 and con- 
tinuum. The HCN and HCO + 3-2 emission is elongated in the 
W/NW-E/SE direction with a clear velocity gradient along its 
major axis. The major axis of the HCO + and HCN 3-2 emission 
is close to the axis of the continuum emission - and perpendic- 
ular to the directio n toward the Herbig-Haro objec ts (Fig . IB. 1 1 
Grosso et alj|200l|) and the thermal continuum jet (Girart et all 
20001) . 



B.2. TMR1 andTMCIA 

Images of TMR1 u sing the HST/NICMOS at 1.6-2.1 p.m by 
Terebev et al. (1998) show s triple system consisting of a close 
binary (TMR-1A and TMR-1B) with a projected separation of 
0.3" and a third source TMR-1C offset from the binary by 
about 10" but connected with this through an extended fila- 
ment of near-infrared nebulositv. lTerebev et al. ( 1998, 2000) dis- 
cuss the option that TMR-1C is a runaway giant planet, but 
conclude based on near-infrared spectroscopy that it is most 
likely is a background star. TMC1A is an emb edded protostar 
with an associated scattering n ebulosity (e.g., Chandler et all 
1996; iHogerheiide et all 1 19981) aligned with blue-shifted out- 
flow emission. 

Both TMR1 and TMC1A are detected in continuum, to- 
gether with HCO + and HCN 3-2. Both show clear evidence 
of the line emission being offset from the continuum peaks 
toward the outflow nebulosity observed in the near-infrared 
HST/NICMOS images - similar to the case for GSS30-IRS1. 
The offsets in the HCO + 3-2 emission towa rd the two sources 
are sim ilar to those observed in HCO + 1-0 by Hogerheiide et al. 
(fl998T) . 



B.3. GSS30 

The GSS30 region is the most complex of the studied fields asso- 
ciated with a clear near-infrared nebulosity. At least three young 
stellar objects (GSS30-IRS1, -IRS2 and -IRS3) are seen at in- 
frared wa velengths, located wi thin the SMA primary beam field 
of view, ( Weintraub et all 1 19931) and observed by Spitzer (e.g., 
lJ0rgensen et all 120081) . GSS30-IRS 1 is located at the center of 
this outflow cone and has a mid-infrared SED slope consistent 
with a Class I protostar. GSS 30-IRS3 (LFAM1) was det ected 
in hi gh-resolution image s at 2.7 mm (Zhang et al., 1997) and 
6 cm dLeous et all 1 19911) . GSS30-IRS2 is likely a more evolved 
(T Tauri) young stellar object also detected at 6 cm. 

The SMA continuum observations identify one protostar 
GSS30-IRS3 (LFAM1) whereas the two other protostars remain 
undetected down to the noise level of « 3 mJy beam 1 . The 
sa me was the case fo r the similar in the 2.7 mm observations 
by|Zhang et al] d 1997b . The single-dish continuum emission to- 
ward GSS30 is of ten attributed solely to GSS30-IRS1 (e.g., 
IZhang et all 119971) . but appears to be peaking at GSS30-IRS3 
consistent with the suggestion that this source is more deeply 
embedded. 

The HCO + line emission shows a distinct peak at the loca- 
tion of GSS30-IRS1 with a narrow line detected toward GSS30- 
IRS2 (Fig.©. HCN is also clearly detected toward GSS30, again 
toward GSS30-IRS1, but offset from the HCO + 3-2 emission. 
Toward IRS 1 , both HCO + and HCN appear to trace one side of 
the outflow cone directed toward us. 



B.4. WL 12 

WL 12 is a "standard" Class I YSO: it is neither a binary, nor 
does it have a spectacular outflow. It is clearly detected in the 
interferometric continuum data, but not in the HCO + or HCN 3— 
2 maps. It is identified as a core i n SCUBA maps and inc luded 
in the list of embedded YSOs by |j0rgensen et all ([2008 ) with 
red mid-infrared colors chara cteristic of those sources, as well 
as bright HCO + 4-3 emission dvan Kempen etalll2009h . 



J. K. J0rgensen et al.: The mass evolution of embedded protostars, Online Material p 4 
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Fig. B.l. Upper panel: The region of the IRS 54 outflow as well as the near-infrared Herbig-Haro objects from Spitzer observations 
with 3.6 pm in blue, 4.5 pm in green and 8.0 pm in red. In addition we show the suggested direction of the precessing IRS 54 
outflow (dotted line) and the direction to IRS 43 and IRS 44 (dashed arrow). Lower panels: images of the near-i nfrared Herbig-Haro 
knots from 3.6-8.0 pm (IRAC1-4) and 24 pm (MIPS1). The location of the Herbig-Haro knots suggested bv lGrosso et alj OOOlT) 
to be associated with the IRS 54 outflow are shown with the crosses, whereas the knots proposed to be associated with the IRS 43 
outflow are shown with plus-signs in the first panel. 



B.5. IRS 54 consistent with the suggestion that it is a source in transition be- 

tween the embedded and T Tauri stages with M env ~ Mdi s k- 

IRS 54 is found to be a binary with a separ ation of about 7" in 
deep near-infrared w avelength observations (lHaisch et all 120041 
Duc hene et al.L 120041) . possibly also responsible for the preces- 
sion of the outflow observed in the mid-infrared Spitzer images 
(Fig. IB. 11 1. IRS 54 is the only source in our sample, which is 
not detected in either continuum and only shows a tentative de- 
tection in HCO + 3-2. It is associated with a faint SCUBA core 
(peak of 0.1 Jy beam -1 ; total integrated flux in a 40" radius 
aperture of 0.4 Jy) be low the threshold for core detection in 
J0rg ensen et al] d2008l) . It is also not identified as a separ ate mil- 
limeter core in the Bolocam maps of Young et al. (2006). Based 
on its HCO + 4-3 line emission, Ivan Kempen et al .1 (|2009) clas- 
sify it as a "late Stage 1" source in transition to the T Tauri stage. 



B.6. IRS 63 and Bias 29 

The r esults for IRS 63 and Elias 29 are discussed in more de- 
tail in Lomme n et al] (12008). In summary, the two sources are 
detected in both continuum and HCO + 3-2. IRS 63 shows the 
strongest continuum emission by about a factor 8 whereas the 
Elias 29 is stronger in HCO + 3-2 also by about a factor 8. Where 
Elias 29 clearly is associated with strong extended continuum 
emission, witnessed both by its brightness profile in the SMA 
data and by comparison to single-dish observations, IRS 63 is 
dominated by emission from a compact, unresolved component, 



